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INTRODUCTION 



This manual contains descriptions of the programs available in the Time-Sharing Library. 
Each program is described using a standard format - referred to as a program abstract. 



Program Abstract 

This specifies the purpose of the program, gives instructions for its use, and demonstrates 
both of these via a sample problem. This is fully stated with the necessary input and 
output provided as it would appear to a user at a terminal. 



Locating a Program Abstract 

This is a simple procedure. Programs are grouped into broad classes: 
B = Business and Financial 
E = Engineering 

M = Mathematical and Scientific 
S = Statistical 
T = Tutorial 

Within each class the program is given a number in the range 001 - 999 This number 
simply indicates the point in time when the program was entered in the library,, Hence 
the latest program will have the highest number in its class. 

Within the manual, the classes are in separate sections; each section has its class letter 
printed at the edge of the page. Hence to find a particular program, it is necessary only 
to go to the relevant section - business, engineering, etc. - and locate the program abstract 
in its numerical sequence. 

The class and number of a program are called the program code. This is printed at the 
top of each page for easy reference. 

For example: The program ANNUIT*** has the program code BOOl. To locate the abstract, 
turn to section (B) and then the first program. 
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In practice, programs aie r( fened to by their names 'atlier liian by program coiie. To 
simplify locating a program an index is provided listing both program names and prof rain 
cudes„ This is in order by program code, reflecting the chror )lo;-ical sequence of progr, in 
entries,, This also makes updating' (ho index very simple. 

P iOGRAM INDEX: 

Note that all program names :ire assumed to have three asterisks following) 



Program 
Code 



Name 



BUSINESS AND FINANCIAL: 



Description 



B001 
B002 
B003 



ANNUIT 
TRUINT 
SAVING 



Computes annuities, loans, mortgages 
Computes true annua] interest on installment loan 
Computes value of saving plan a( Compound 
interest 



ENGINEERING: 



E001 


LPFILT 


E002 


LCVSIC 


E003 


GCVSIZ 


E004 


SCVSIZ 


EDO 5 


BEMDES 



To design lowpass filters 
To compute liquid value coefficients 
To compute gas/vapor valve coefficients 
To compute steam valve coefficients 
Steel beam design program 



MATHEMATICAL AND SCIENTIFIC: 
M001 ZEROES 



M002 


CIRCLE 


M003 


SIMEQN 


M004 


ROOTER 


M005 


NUMINT 


M006 


SPHERE 


M007 


TRIANG 


M008 


CURFIT 


M009 


XYPLOT 


MO 10 


PLOTTO 


M011 


TvVOPLO 


MO 12 


LINPRO 


M013 


TMFCEV 


M014 


MULFIT 


M015 


CPMONE 



To find zeroes of a function and relative maxima 

and minima 
To divide a circle into (n) given parts 
To solve (n) linear equations in (n) unknowns 
To find roots of polynomial (Bairstow's method) 
To evaluate definite integral (Gauss' method) 

To solve spherical triangles 

To solve for major features of a triangle 

To determine best fit of G curves to data 

To plot single function of X 

Plots 1 to 6 functions simultaneously 

Plots two functions oi X simultaneously 
Maximize.; objective function (2-phase method) 
Evaluates time functions (exp. + sin/cos) 
Supplies statistical analysis of up to 61 observations 

(1-6 variables) 
CPM input edit program (see CPMONX***) 
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STATISTICAL: 




S001 


MANDSD 


S002 


UNISTA 


S003 


STATAN 


S004 


PROVAR 


S005 


LINREG 


S006 


COLINR 


S007 


CONBIN 


S008 


CONDIF 


S009 


CONLIM 


S010 


BICONF 


SOU 


BITEST 


SO 12 


BINDIS 


S013 


SEVPRO 


SO 14 


ONEWAY 


S015 


TESTUD 



(List program STADES*** for general info on 
several of programs below) 

Computes mean, variance, std. deviation for data 
Description of univariate data 
Gives 34 measures for univariate data 
Calculations for normal and Student's distribution 
Linear regression: three curve types 

Confidence limits on linear regression 
Confidence limits on proportion using normal 

approximation 
Confidence limits on difference between 2 means 
Confidence limits on sample mean value 
Confidence limits on proportion using binomial 

Test of binomial proportion 
Computes binomial probability distributions 
Applies chi-square test to sample proportions 
A variance analysis with equal sample sizes 
Tests population means using sample statistics 



Note: 

The above programs have all been reviewed and tested 
assumed. 



however, no responsibility can be 
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BUSINESS & FINANCIAL 
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Name ANNUIT * * * 
Code B001 



PROGRAM ABSTRACT 



PURPOSE 

To perform the calculations necessary for determining both payment and withdrawal annuities. 

INSTRUCTIONS 

To use this program supply values for the variables as required by the problem. 

Variables are denoted as follows: 

N = Number of periods 

P = Original principal amount 

A = Total amount at end of n periods 

I = Interest rate per period, in percent 

R = Amount of payment/withdrawal each period 

For a payment annuity, you may give any three of N, A, I, R and find the fourth. 

For a withdrawal annuity, you give any three of N, P, I, R and find the fourth. 

For loan or mortgage, use the withdrawal annuity option. 

Alter each case, you may choose one of the following alternatives: 

1 = Another case, same unknown variable 

2 = Another case, different unknown 

3 = Another case, other type of annuity 

4 = Total interest paid over n periods 

5 = Table of withdrawals, principal, and interest 

6 = Stop the program 

Additional instructions maybe found in the listing. 

SAMPLE PROBLEM 

Payment Option: 

Assume you intend to deposit $35 a monlh in the local bank and want to know how much you 
would have after four years. The bank pays 4-7/8% a year and compounds quarterly. 
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ANN U IT * * * 

IJ001 

I'aee 2 



Fur this program you must convert everything to a quarterly basis. 

N = 4 . 4 = 16 pel iuds 

A = varialile oi iulercsl 

1 = 4.875,4 = 1.2l875/> i lod 

R = 35 . 4 = lO'i/l't-riod 

The al)o\ e values of N, 1 and R will be entered during' the runnine, of the program. 



NOTE: 



You must supply the portions of the printout which are underlined. 

The answer ($1842.65) is printed out after you select Ihe proper options at lines 2 and 
3 and enter (in proper sequence) N, I and R at line 4. 

The answer does not account for any simple interest that might have been paid on $35 
deposits prior to the first compounding period. 

The total interest (4) and table (5) alternatives do not apply to the payment option. 



Withdrawal Option: 

Determines the monthly principal and interest payments on a $20,900 mortgage for 30 years 
at 5-1/2% per year. 

To use this program you must convert everything to a monthly basis: 

N = 30 . 12 = 360 periods 

P= 20900 

1= 5.5/12= .458333333/period 

R = variable of interest 

The above values of N, P, and I, will be entered during the running of the program. 

The complete run will be as shown in the sample solution. 

1. The answer ($118. 668) is printed out when the proper options at lines 2 and 3 are 
entered (in proper sequence) N, P, and I, at line 4. 

2. The total interest and a payments table was obtained by answering 4 and 5 respec- 
tively to the "another case" question. The "S" key was used to stop the table 
print-out at period 13. 
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ANNUIT * * * 

B001 

Page 3 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 
RUN 

ANNUIT 

ANNUITY CALCULATI0N 

WANT INSTRUCTIONS FIRST tUYES, 2=N0]?2 
WHICH ANNUITY TYPE INPAYMENT, 2 = WITHDRAWAL]71_ 
WHICH VARIABLE IS UNKN0WN [I =N,2=A,3=I, 4=R ]?2 
WHAT ARE N[INTEGER1,I(PCT),R(S17 l6 t 1. 21875, 105.00 
AM0UNT AT END 0F N PERI0DS = A = 1842.65 
AN0THER CASE (TYPE C0DE NUMBER )?6 
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AMNUIT * * * 

B001 

Page 4 



RUN 



ANNUIT 



ANNUITY CALCULATION 

WANT INSTRUCTIONS FIRST llrYES, 2=N0]?2_ 
WHICH ANNUITY TYPE INPAYMENT, 2=WITHDRAWAL )? 2 
WHICH TRIABLE IS UNKNOWN tl = N,2=P,3=I, 4=R ]? 4_ 
WHAT ARE NCINTEGER),P[$],IfPCT]? 360, 20900, .458333333 
WITHDRAWAL EACH PERIOD : R : 118.668 
ANOTHER CASE [TYPE CODE NUMBER J7 4 




T0TAL INTEREST PAID r 21820.5 

ANOTHER CASE [TYPE CODE NUMBER ]75 



PERIOD 



I 

2 

3 

4 

5 

6 

7 

3 

9 

10 

11 

12 

13 

14 

15 

16 



20900 



'RINCIPAL 


INTEREST 


PRINC BAL 
20900 


INT TO DATE 


22.88 


95.7879 


20877.1 


95.7879 


22.98 


95.6879 


20854.1 


191.476 


23.09 


95.5779 


20831 .1 


287.054 


23.19 


95.4779 


20807.9 


382.532 


23.3 


95.3679 


20784.6 


477.9 


23.41 


95.2579 


20761.2 


573.157 


23.51 


95.1579 


20737.6 


668.315 


23.62 


95.0479 


20714. 


763.363 


23.73 


94.9379 


20690.3 


858.301 


23.84 


94.8279 


20666.5 


953.129 


23.95 


94.7179 


20642.5 


1047.85 


24.06 


94.6079 


20618.4 


1142.45 


24.17 


94.4979 


20594.3 


1236.95 


24.28 


94.3879 


20570. 


1331.34 


24.39 

•2T 


94.2779 


20545.6 


1425.62 
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Name TRUINT * * * 
Code B002 



PROGRAM ABSTRACT 



PURPOSE 

To calculate the true annual interest rate charged on an instalment loan. 

INSTRUCTIONS 

To use this program simply supply values for the four variables: 
A = AMOUNT OF LOAN [IN $] . 
P = AMOUNT OF EACH PAYMENT [$] . 
N = THE TOTAL NUMBER OF PAYMENTS DUE. 
K = THE NUMBER OF PAYMENTS DUE IN ONE YEAR. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Determine the true interest charged on these two loans: 

Loan 1 is for $600 for 21 months with a payback of $31. 99 per month. Loan 2 is for $1400 
for 30 months with a payback of $54.90 per month. The bank's stated rate is 6%. The input 
data would be: 



Loan 1 



A = Amount of Loan = $600. 00 

P = Amount of each payment = $31. 99 

N = Total Number of Payments =21 

K = Number of Payments in one year = 12 



Loan 2 



A = $1400.00 
P = $54.90 
N = 30 
K = 12 
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TRUINT * * * 

B002 

Patre 2 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

Note: 

All user supplied information is underlined 

TRUINT 

THIS PROGRAM WILL CALCULATE THE TRUE ANNUAL INTEREST 
RATE CHARGED 0N AN INSTALLMENT L0AN. Y0U SUPPLY THE 
VALUES 0F F0UR VARIABLES: A = AM0UNT 0F L0AN UN $], 
P : AM0UNT 0F EACH PAYMENT I$J, N = THE T0TAL NUMBER 
0F PAYMENTS DUE, AND K = THE NUMBER 0F PAYMENTS DUE 
IN 0NE YEAR. WHAT ARE A, P, N,K? 600.00, 31.99, 21, 12 

THE TRUE ANNUAL INTEREST RATE - 12.61 

AN0THER CASE? tTYPE 'S' T0 ST0P N0W1. 
WHAT ARE A,P,N,K? 1400.00, 54.90, 30, 12 

THE TRUE ANNUAL INTEREST RATE : 12.98 

AN0THER CASE? tTYPE 'S' T0 ST0P N0W). 
WHAT ARE A,P,N,K?S 
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Name SAVING * * * 
Code B003 



PROGRAM ABSTRACT 



PURPOSE 

To calculate the amount of money that would accumulate after N years at an annual interest 
rate R compounded T times per year, when the initial amount is P and an amount D is 
added at the beginning of each subsequent year. 



INSTRUCTIONS 

To use this program enter the data as requested. 

NOTE: 

P and D must be given in dollars, N and T in integers, and R as a percentage. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Several cases are used based on the input data underlined in the sample solution. 
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SAVING * * * 
B003 
Pase 2 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 
RUN 
SAVING 

THIS PR0GRAM CALCULATES THE AMOUNT 0F M0NEY THAT 
W0ULD ACCUMULATE AFTER N YEARS AT AN ANNUAL INTEREST 
RATE R C0MP0UNDED T TIMES PER YEAR, WHEN THE INITIAL 
AM0UNT IS P AND AN AM0UNT D IS ADDED AT THE BEGINNING 
0F EACH SUBSEQUENT YEAR. N0TE THAT P AND D ARE GIVEN 
IN DOLLARS, N AND T MUST BE INTEGERS, AND R IS GIVEN 
AS A PERCENTAGE. 

WHAT ARE P,D,N,T,R? 6000.00.200.00,2Q,4,4 

AFTER 20 YEARS, 6000 D0LLARS INVESTED AT 4 

PERCENT C0MP0UNDED 4 TIMES PER YEAR, WITH THE 

ADDITI0N 0F 200 D0LLARS PER YEAR, YIELDS A T0T4L 
0F 19093.4 D0LLARS. 

TYPE 'S* T0 ST0P N0W, 0R ELSE 

WHAT ARE P, D,N t T,R? 6000.00,0, 20, 4, 4 

AFTER 20 YEARS, 6000 D0LLARS INVESTED AT 4 

PERCENT C0MP0UNDED 4 TIMES PER YEAR, WITH THE 

ADDITI0N 0F D0LLARS PER YEAR, YIELDS A T0TAL 

0F 13300.3 D0LLARS. 

TYPE 'S' T0 ST0P N0W, 0R ELSE 

WHAT ARE P,D,N,T,R? 6000.00,0,20,1,4 

AFTER 20 YEARS, 6000 D0LLARS INVESTED AT 4 

PERCENT C0MP0UNDED 1 TIMES PER YEAR, WITH THE 

ADDITI0N 0F D0LLARS PER YEAR, YIELDS A T0TAL 

0F 13146.7 D0LLARS. 

TYPE 'S' T0 ST0P N0W, 0R ELSE 

WHAT ARE P,D,N,T,R? 7000.00,0,20,4,5.25 

AFTER 20 YEARS, 7000 D0LLARS INVESTED AT 5.25 

PERCENT C0MP0UNDED 4 TIMES PER YEAR, WITH THE 

ADDITI0N 0F D0LLARS PER YEAR, YIELDS A T0TAL 

0F 19867.4 D0LLARS. 

TYPE 'S* T0 ST0P N0W, 0R ELSE 
WHAT ARE P,D,N,T.R? 6000,0,3,4,5.50 

AFTER 3 YEARS, 6000 D0LLARS INVESTED AT 5.5 

PERCENT C0MP0UNDED 4 TIMES PER YEAR, WITH THE 

ADDITI0N 0F D0LLARS PER YEAR, YIELDS A T0TAL 

0F 7068.41 D0LLARS. 

TYPE *S' T0 ST0P N0W, 0R ELSE 
ilHAT ARE P,D,N,T,R7S_ 
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ENGINEERING 




-GE-400 SERIES TIME-SHARING 



Name LPFILT * * * 
Code E001 



PROGRAM ABSTRACT 



PURPOSE 



To design low pass filters using constant K prototype T section and M derived |. 1=0. 6j 
termination L sections. Up to nine additional M derived T sections may be included to give 
high attenuation at specified frequencies in the stop band. 



INSTRUCTIONS 

To use this program, enter data as: 

10 DATA R 7 C, N, F[l], F [2] ,..., F£n] 

WHERE R = DESIRED CHARACTERISTIC IMPEDANCE IN OHMS 

C = DESIRED CUTOFF FREQUENCY IN CYCLES/SECOND 
N = NUMBER OF ATTENUATORS DESIRED IN STOP BAND 
F [l] = FREQUENCY FOR ATTENUATOR I 

Then Type RUN 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Design a low pass filter with a 50 ohm impedance. 20 kc cut-off frequency. 2 attenuators 
in stop band, 455KC and 91 KC attenuators in the filter. 
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LPFILT * * * 

E001 

Page 2 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 
NOTE: 

All user supplied information is underlined. 



10 DATA 50. 2E4.2, 455000. 91E3 

sire — ■ 

WAIT. 



LPFILT 



DESIGN F0R DESIRED L0W PASS FILTERj 



0<- 
I 

+-• 
I 

> 

I 
+-• 

I 

1 

+-• 

I 

> 

I 

+-■ 

I 

> 

I 
+- 

I 
0< 



50 0HM LINE 

.424413 MH + .095493 

.63662 MH 

----- .31831 MFD — 

.79539 MH 

3.84758 E-4 MH + .318002 

.735662 MH 

9.85050 E-3 MH + .310527 

.626891 MH 

.424413 MH + .095493 

50 0HM LINE 



•>0 

I 



MFD 



MFD 



MFD 



MFD 



■-♦ 

I 
•>0 



TERMINATING SECTIONS GIVE MAXIMUM ATTENUATI0N AT 25000 
CPS IN ADDITI0N T0 THE SPECIFIED ATTENUAT0R FREQUENCIES. 
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Name LCVSIC * * * 
Code E002 



PROGRAM ABSTRACT 



PURPOSE 

To calculate liquid valve coefficients and the required valve rangeability. 

INSTRUCTIONS 

Enter input data using the following format : 

TYPE: 20 DATA [ENTER DATA IN ORDER. . .MIN. PRESS. DROP, 
MAX. PRESS. DROP, MIN. FLOW AND MAX. FLOW. 

TYPE: 25 DATa[eNTER ONE OR MORE VALUES OF FLUID SPECIFIC 
GRAVITY, FOR EXAMPLE, EXPECTED HIGH AND LOW. 

Then type RUN 

NOTE: 

To size more valves, enter new data in lines 20 and 25, then type RUN. 

Reference is the voluntary standard (Nov. , 1961) of the Fluid Controls Institute, Inc. 
No viscosity corrections are available yet. Consult valve manufacturers about viscous or 
non-newtonian fluids. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Determine the valve coefficients which satisfy the data used in the sample solution. 
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LCVSIC * * * 
E002 
Patre 2 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 



20 DATA 30, 70, 3.5, 12.5 

25 DATA .9. "1. 1.1 

RUN 

WAIT. 



LCVSIC 1U40 

***THE IDENTITY 0F THE VALVE IS 
WHEN SPECIFIC GRAVITY = .9 

FL0W.GPM PRESS. DR0P.PSI VALVE C0EFF. 

3.5 30 .606218 

12.5 30 2.16506 

3.5 70 .396863 

12.5 70 I. 41737 

RATI0 0F HIGH/L0W C0EFFICIENTS: VALVE RANGEABILITY 

WHEN SPECIFIC GRAVITY = I 

FL0W.GPM PRESS. DR0P.PSI VALVE C0EFF. 

3.5 30 .63901 

12.5 30 2.28218 

3.5 70 .41833 

12.5 70 1.49404 

RATI0 0F HIGH/L0W C0EFFICIENTS: VALVE RANGEABILITY 

WHEN SPECIFIC GRAVITY =1.1 

FL0W.GPM PRESS. DR0P.PSI VALVE C0EFF. 

3.5 30 .670199 

12.5 30 2.39357 

3.5 70 .438748 

12.5 70 1.56696 

RATI0 0F HIGH/L0W C0EFFICIENTS: VALVE RANGEABILITY 

0UT 0F DATA IN 40 
TIME: 1 SECS. 
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Name GCVSIZ * * * 
Code E003 



PROGRAM ABSTRACT 



PURPOSE 

To calculate gas and vapor (but not steam) control valve coefficients and the required valve 
rangeability. 



INSTRUCTIONS 

Enter input data using the following format: 

TYPE: 20 DATA |_ENTER IN ORDER MIN. TEMP. IN DEGREES 

F, MAX. TEMP. IN DEGREES F, MOLECULAR WEIGHT. 

TYPE: 25 DATA [ENTER DATA IN GROUPS OF 3, AS MANY GROUPS 

AS YOU NEED, IN THIS ORDER PPH FLOW, INLET 

PSIG, OUTLET PSIG, 2ND FLOW, 2ND PRESS., ETC... 

Then type RUN 

NOTE: 

When pressure drop is critical, the program assumes outlet pressure is half the 
absolute inlet pressure. The program also converts pressure and temperature to 
absolute. 

Reference is equation 4 of the voluntary standard agreed to in Nov. , 1961 by the 
Fluid Controls Institute, Inc. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Determine the control valve coefficients which satisfy the data used in the sample solution. 
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SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 



NOTE: 



All user supplied information is underlined. 



20 DATA 470.870,1 12.5 

25 DATA 1350. 145, 80,8000, 145. 80, 1350, 195, 1 05,8000, 195, 105 

RUN 

WAIT. 



GCVSIZ 11:34 

RATI0 0F HIGH/L0W C0EFF.=REQUIRED VALVE RANGEABILITY. 

*** THE IDENTITY 0F THE VALVE IS 



WHEN FL0W IN PPH = 1350 

TEMP IN DEG F INLET PSIG 
470 145 

870 145 



FL0W IN SCFH = 4550.53 

0UTLET PSIG PRESS DR0P PSI VALVE C0EFF. 

80 65. 2.20719 

80 S5. 2.63951 



WHEN FL0W IN PPH = 8000 

TEMP IN DEG F INLET PSIG 
470 145 

870 145 



FL0W IN SCFH : 26966.1 

0UTLET PSIG PRESS DR0P PSI VALVE C0EFF. 

80 65. 13.0796 

80 65. 15.6416 



WHEN FL0W IN PPH = 1350 

TEMP IN DEG F INLET PSIG 
470 195 

870 195 



FL0W IN SCFH = 4550.53 

0UTLET PSIG PRESS DR0P PSI VALVE C0EFF. 
105 90. 1.64843 

105 90. 1.97131 



WHEN FL0W IN PPH = 8000 

TEMP IN DEG F INLET PSIG 
470 195 

870 195 

0UT 0F DATA IN 45 



FL0W IN SCFH = 26966.1 

0UTLET PSIG PRESS DR0P PSI VALVE C0EFF. 

105 90. 9.76849 

105 90. 11.6819 



TIME: 



1 SECS. 
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Name SCVSIZ * * * 
Code E004 



PROGRAM ABSTRACT 



PURPOSE 

To calculate steam control valve coefficients and the required valve rangeability. 

INSTRUCTIONS 

Enter input data using the following format: 

TYPE: 20 DATA [eNTER DEGREES F OF SUPERHEAT IN THE STEAM 
OR ENTER ZERO IF IT IS ZERO. 

TYPE: 25 DATA [ENTER DATA IN GROUPS OF 3 , AS MANY GROUPS 

AS YOU WISH QR NEED, IN THIS ORDER. . .FLOW IN PPH, 
INLET PSIG, OUTLET PSIG ; THEN SECOND FLOW, SECOND 
INLET PRESSURE , SECOND OUTLET PRESS . , ETC 

Then type RUN 



NOTE: 



When pressure drop is critical, the program assumes outlet pressure is half the 
absolute inlet pressure. The program also converts pressures to absolute. 



Reference is equation 8 of the November, 1961 Fluid Controls Institute, Inc. 
voluntary standard. 

To size more valves, enter new data in 20 and 25 and type RUN. If for saturated 
steam, zero in 20 need be entered only once. It will stay in. 



Additional instructions may be found in the listing. 
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SAMPLE PROBLEM 

Determine the valve coefficients which satisfy the data used in the sample solution. 

SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 



20 DATA 463 

25 DATA 1000,106,104,10000,106,104,1000,66,64,10000,66,64 

RUN 

WAIT. 



SCVSIZ 11:17 

RATI0 0F HIGH/L0W C0EFF.:REQUIRED VALVE RANGEABILITY 

*** THE IDENTITY 0F THIS VALVE IS 

WHEN SUPERHEAT IN DEGREES F = 463 



INLET PSIG 


0UTLET PSIG 


PR.DR0P.PSI 


FL0W, PPH 


VftLVE C0EFF 


106 


104 


2 


1000 


28.8154 


106 


104 


2 


10000 


288.154 


66 


64 


2 


1000 


35.3136 


66 


64 


2 


10000 


353.136 



0UT 0F DATA IN 60 



TIME: 



1 SECS. 
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Name BEMDES * * * 
Code E005 



PROGRAM ABSTRACT 



PURPOSE 

To recommend the correct steel beam to be used in any of the following 13 load and support 



cases. 



CASE 1 Uniform distributed load on simple beam. 

CASE 2 Uniform distributed load, beam fixed at one end, supported at other. 

CASE 3 Uniform distributed load, beam fixed both ends. 

CASE 4 Uniform distributed load, cantilever beam. 

CASE 5 Single midpoint load on simple beam. 

CASE 6 Single midpoint load, beam fixed at one end, supported at other. 

CASE 7 Single midpoint load, beam fixed both ends. 

CASE 8 Single midpoint load, cantilever beam. 

CASE 9 Two equal loads, placed symetrically on simple beam. 

CASE 10 Uniform load plus midpoint load on simple beam. 

CASE 11 Uniform load plus midpoint load, beam fixed one end, supported at other. 

CASE 12 Uniform load plus midpoint load, fixed both ends. 

CASE 13 Uniform load plus midpoint load, cantilever beam. 



INSTRUCTIONS 

To use this program type RUN, then respond to the following questions: 

Which load and support case number properly describes your problem (list lines 1 thru 
16 for definitions)? 

What is the length of the span in feet? 

What is the distributed load in pounds per foot? 

What is the amount of each concentrated load in pounds? 

What are the locations of the concentrated loads in feet from the end 1 ? 

Additional instructions may be found in the listing. 
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SAMPLE PROBLEM 

Determine the type of beam required for each of the 13 cases listed in the following table. 



Case 
No. 



PICTURE 



Feet 
Span 



Di st rib 
Load 



Concen 
Load 



From 
End 



10 
11 
12 
13 



DISTRIB 

miiiiiiuui 

A A 

JlllUIIIIIUli 

r lUUUlim.U fc: 



CONC 



% « 

i FROM t 

[VAi end i A 

F — »— *l 

lllilllllllili 



!ii,)i 



.nilililltti. l 



Jtumluiui 



23 




1000 


23 




100 


23 




1000 


23 




300 


23 






23 






23 






23 






2 3 






23 




6000 


23 




15000 


2 3 




1000 


2 3 




500 


CAN NOT 


BE 


SOLVED 


BY THIS 


PROGRAM 



25000 




50000 




50000 




25000 




10000 


(A)4' 


10000 


(B)9' 


150000 




4 00000 




25000 




10000 
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SAMPLE SOLUTION 

The following sample solutions are copies of the printout that will appear at your teletype- 
writer. 



NOTE: 

All user supplied information is underlined. 
RUN 

BEMDES 

PR0GRAM FOR SELECTION OF STEEL BEAMS. 

WHICH LOAD AND SUPPORT CASE NUMBER PR0PERLY DESCRIBES Y0UR 

PROBLEM (LIST LINES 1 THRU 16 FOR OEFINITI0NS) ? 1 

WHAT IS THE LENGTH 0F THE SPAN IN FEET t gi 

WHAT IS THE DISTRIBUTED L0AD IN P0UNDS PER F00T ? 1000 

USE A 14 WF 30 



RUN 

WAIT. 



BEMDES 

PR0GRAM F0R SELECTION 0F STEEL BEAMS. 

WHICH L0AD AND SUPP0RT CASE NUMBER PR0PERLY DESCRIBES Y0UR 

PR0BLEM CLIST LINES 1 THRU 16 F0R DEFINITIONS) ? £. 

WHAT IS THE LENGTH 0F THE SPAN IN FEET ? 2J 

WHAT IS THE DISTRIBUTED L0AD IN P0UNDS PER F00T ? 100 

USE AS JR 6.5 



RUN 

BEMDES 

PR0GRAN F0R SELECTION 0F STEEL BEAMS. 

WHICH L0AD AND SUPPORT CASE NUMBER PROPERLY DESCRIBES YOUR 

PROBLEM [LIST LINES 1 THRU 16 FOR DEFINITIONS! ? i_ 

WHAT IS THE LENGTH OF THE SPAN IN FEET ? 23 

WHAT IS THE DISTRIBUTED LOAD IN POUNDS PER FOOT ? 1000 

USE A 12 WF 27 
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RUN 
WAIT. 



BEMDES 

PROGRAM FOR SELECTION 0F STEEL BEANS* 

WHICH LOAD AND SUPPORT CASE NUMBER PROPERLY DESCRIBES YOUR 

PROBLEM (LIST LINES I THRU 16 FOR DEFINITIONS] ? 8 

WHAT IS THE LENGTH OF THE SPAN IN FEET ? 23 

WHAT IS THE AMOUNT OF EACH CONCENTRATED LOAD IN POUNDS ? 23000 

USE A 24 WF 84 



RUN 

BEMDES 

PROGRAM FOR SELECTION OF STEEL BEAMS. 

WHICH LOAD AND SUPPORT CASE NUMBER PROPERLY DESCRIBES YOUR 

PROBLEM (LIST LINES I THRU 16 FOR DEFINITIONS] ? S 

WHAT IS THE LENGTH OF THE SPAN IN FEET ? 23 

WHAT IS THE AMOUNT OF EACH CONCENTRATED LOAD IN POUNDS ? 10000 

WHAT ARE THE LOCATIONS OF THE CONCENTRATED LOADS IN FEET 

FROM THE END ? 4 9 

USE A 30 WF 108 



RUN 

BEMDES 

PROGRAM FOR SELECTION OF STEEL BEAMS. 

WHICH LOAD AND SUPPORT CASE NUMBER PROPERLY DESCRIBES YOUR 

PhOBLEM (LIST LINES 1 THRU 16 FOR DEFINITIONS] ? 13 

WHAT IS THE LENGTH OF THE SPAN IN FEET ? 23 — 

WHAT IS THE AMOUNT OF EACH CONCENTRATED LOAD IN POUNDS ? 10000 

WHAT IS THE DISTRIBUTED LOAD IN POUNDS PER FOOT ? 500 

USE A 24 WF 76 



NOTE: 



Solutions for sample problem cases 1, 2, 3, 8, 9, and 13 provide at least one example 
of each unique data situation. 
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PROGRAM ABSTRACT 



PURPOSE 

To locate "interesting" values of X for any function of X. Specifically to local ■ the values 
of X at which relative maximums and minimums of F [x] occur, and the values of X for 
which F [x] is zero [i.e. , the zeroes or roots of the function] . 



INSTRUCTIONS 

To use this program, enter data using the following format: 

100 DATA <XMIN>,<XMAX>,<ACC>,<INCR> 

200 LET Y «= <A FUNCTION OF X> 

RUN 

NOTE: 

XMIN and XMAX define the interval in which interesting values of X are to be sought, 
ACC is the accuracy [in number of significant figures] to which the ZEROES of X 
and the maximum and minimum values of F [x] are to be estimated, INCR is the 
number of increments into which the total interval is to be divided for search purposes 
[try 50 to start] , and the function of X is any legitimate basic-language expression 
involving the variable X. 

Additional instructions may be found in the listing. 



SAMPLE PROBLEM 

To illustrate how this program is used consider the sine function, it crosses the axis at 
degrees and every 180 degrees thereafter, has maximums at 90 degrees and every 180 
degrees thereafter and has minimums at 270 degrees and every 180 degrees thereafter. 

To use this program to arrive at the same conclusion enter the following data: 

100 DATA 0, 725, 3, 200 

200 LET Y = SIN(X*3. 14159265/180) 
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NOTE: 



Since BASIC assumes all trig functions are in radians, you must convert degrees to 
radians by multiplying X by Pi over 180. 

The slight errors at the maximums and minimums are due to the binary nature of the 
machine- -Pi could not be represented exactly and it was stepping through the function 
in fractional steps. 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 



100 DATA 0, 725, 3, 200 

200 LET Y ; SI N(X*3 . 1 4 159265/ 180 ) 

RUN 



ZER3ES 



P3INT-TYPE 


F 


ZER3 


0. 


MAX 


1. 


ZER0 


0. 


!*!IN 


-1 . 


ZER0 


0. 


MAX 


1. 


ZER0 


0. 


MIN 


-1. 


ZER3 


0. 



FtXl 





90.002 

180.004 

270.006 

360.008 

450.01 

540.012 

630.014 

719.987 
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PROGRAM ABSTRACT 



PURPOSE 

To divide circles into any number of equal parts, giving the angles in Decimal Degrees and 
Degrees, Minutes, and Seconds, and to calculate the Horizontal and Vertical Distances from 
the center to the point on the circumference. 



INSTRUCTIONS 

To use this program enter the input data in the following format: 
10 DATA N, R, 



NOTE: 

N = the number of parts 
R = the radius of circle 

As many cases as desired may be entered in the same way by continuing the data list. 

After all the data is entered type RUN. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Divide 3 sheet metal disks of 18. 351 inches radius into equal seqments of 3, 7, and 15 
parts. To use this program enter the following data: 

10 DATA 3, 18.351, 7, 18.351, 15, 18.351 
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SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

READY. 

10 DATA 3, 18.331, 7, 18.351, I?, 16.351 

RUN 
WAIT. 

CIRCLE 



DIVISION OF 



CIRCLE 



CASE NUMBER t 1 

ANGLE.... 

INDEX DECIMAL DEQ MIN 



1 120 

2 240 

3 360 



number of parts: 3 
""sec 



120 
240 
360 



RADIUS: 18.351 

CO-ORDINATES 

HORIZ VERT 



■9.1755 

•9.17551 

18.351 



15.8924 
-15.8924 




CAS 


:E NUMBER: 


2 


NUMBER OF PARTS: 7 


RADIUS: 


18 


.351 












. . .. .CO-ORDt m*tit«. 


INDEX 


DECIMAL 


DEQ 


MIN 


SEC 


• • 9 * ♦ w \J Vi\ W f 

HORIZ 




VERT 


1 


51.4286 


51 


25 


42.9 


11.4417 




14.3474 


2 


102.857 


102 


51 


25.7 


-4.08348 




17.8909 


3 


154.286 


154 


!7 


8.6 


-16.5337 




7.96221 


4 


205.714 


205 


42 


51.4 


-16.5337 




-7.96219 


5 


257.143 


257 


8 


34.3 


-4.08349 




-17.8909 


6 


308.571 


308 


34 


17.1 


11.4417 




-14.3474 


7 


360. 


360 








18.351 
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CASE NUMBER t S 



NUMBER OF PARTS: 15 





DECIMAL 


DE6 






NDEX 


MIN 


SEC 


1 


24 


24 








2 


48 


48 








3 


72 


72 








4 


96 


96 








5 


120 


120 








6 


144 


144 








7 


168 


168 








8 


192 


192 








9 


216 


216 








10 


240 


240 








11 


264 


264 








12 


288 


288 








13 


312 


312 








14 


336 


336 








15 


360 


360 









RADIUS* 18.351 

CO-ORDINATES 

HORIZ VERT 



16.7645 


7.46402 


12.2792 


13.6374 


5.67077 


17.4528 


1.9182 


18.2505 


9.1755 


15.8924 


14.8463 


10.7865 


17.95 


3.8154 


17.95 


-3.81538 


14.8463 


•10.7864 


9.17551 


-M.8924 


1.91821 


-I ,.2505 


5.67076 


-17.4528 


12.2792 


-13.6375 


16.7645 


•7.46404 


18.351 
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PROGRAM ABSTRACT 



PURPOSE 

To solve systems consisting of N linear equations in N unknowns. 

INSTRUCTIONS 

To use, this program, enter coefficients in the equations in data statements, starting in 
statement number 1 with the first coefficient of the first equation, and ending with the N-TH 
coefficient of the N-TH equation. All zero coefficients must be entered in their proper 
place. The right-side constant terms of the equations are then entered in subsequent data 
statements. If additional cases with the same coefficient matrix but different right sides 
are to be run, they may all be run at once by simply entering additional data statements 
with the right-side values of the other cases. A zeroth data statement [preceding all the 
above] is used to specify the number of systems to be solved, and the number of equations 
[and hence variables] in the system. Thus, the two systems: 



3X 


4 


5Y 


- 


2Z 


= 


9 


7X 


+ 


Y 






= 


-3 


X 


- 


7Y 


+ 


9Z 


= 


14 



3X 


+ 


5Y 


- 


2Z 


= 


19 


7X 


+ 


Y 






= 


-3 


X 


- 


7Y 


+ 


9Z 


= 


8 



Could be solved by typing: 

DATA 2,3 

1 DATA 3,5,-2,7,1,0,1,-7,9 

2 DATA 9,-3,14,19,-3,8 

RUN 

Solutions are proofed by multiplying the vector by the original coefficient matrix. 
Additional instructions may be found in the listing. 
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SAMPLE PROBLEM 

Solve the following sets of simultaneous equations: 



25W 


-2.7X 


+38. 7Y 


+ Z = 115.76 


10 




18X 


- 5.8Y 


- 3Z = 11.5 


20 


6W 


+ 9X 


+ 18Y 


= 55 


30 


11W 




+ 87Y 


+41Z = 17 


40 



The extreme right-hand column of numbers is a second set of right-hand values. 



To use this program to solve the two sets, prepare the following data tape: 

DATA 2, 4 

1 DATA 25, -2.7, 38.7, 1, 0, 18, -5.8, -3 

2 DATA 6, 9, 18, 0, 11, 0, 87, 41 

3 DATA 115.76, 11.5, 55, 17, 10, 20, 30, 40 



In theO line, the 2 means two sets of equations and the 4 means each set consists of four 
equations. In the 1 and 2 lines, the 0's are necessary to represent any missing (zero coef- 
ficient) terms in the set. Line 4 contains the right-hand values for both sets. 



NOTE: 



Supply all portions of the print-out which are underlined. 



The answers for the first set would be read thus: 



W = 1.25315 
X = .578122 



Y = 2.34878 
Z = -4.90557 



The proof lines merely use the calculated values to evaluate the left-hand side of 
each equation to prove that they give the right-hand value. 




GE-400 SERIES TIME-SHARING 



SIMEQN * * * 
M003 
Page 3 



SAMPLE SOLUTIONS 

This sample solution is a copy of the printout that will appear at your teletypewriter. 



NOTE: 



All user supplied information is underlined. 



DATA 2. 4 

1 DATA 25~-2.7. 38.7. I. 0. 18. -3.8. -3 

2 DATA 6. 9. 18. 0. 11, 0. 87. 41 

3 DATA 115.76, 11.5, 55, 17, lO, 20, 30, 40 
ffUfi 



SIMEQN 

S0LUTI0N F0R LINEAR SYSTEM 0F ORDER 4 

INDEXt 
1 2 3 



S0LUTI0N VECT0R F0R CASE I 
1.25315 .578122 



2.34878 



-4.90557 



PR00F 0F S0LUTI0N F0R CASE 1 
115.76 11.5 55. 



17. 



S0LUTI0N VECT0R F0R CASE 2 
-1.95156 1.3087 



1.66283 



-2.02925 



PR00F 0F S0LUTI0N F0R CASE 2 
10. 20. 30. 



40. 
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PROGRAM ABSTRACT 



PURPOSE 

To find the roots of Polynomials using Bairstow's method. 

INSTRUCTIONS 

To use this program, enter data in data statements as follows: 

10 DATA N, A(N) , A(N-l) , A(l) , A(0) 

WHERE N IS THE ORDER OF THE POLYNOMIAL AND A(I) IS THE COEFFICIENT OF THE 
I-TH DEGREE TERM. MORE THAN ONE DATA LINE MAY BE USED TO SUPPLY 
COEFFICIENTS FOR ONE POLYNOMIAL, AND ADDITIONAL POLYNOMIALS MAY BE 
SOLVED ON A SINGLE RUN BY SUPPLYING DATA FOR THEM ON SUBSEQUENT DATA 
LINES (NOT BEYOND LINE 299). 

NOTE: 

There are a few types of polynomials which this program is unable to solve. Program 
will so indicate and go on to next case. Also, for high order polynomials, program 
may require many iterations, so running time can be high. 

Additional instructions may be found in the listing. 



SAMPLE PROBLEM 

Determine the roots for the following 6 polynomials: 



6 + 3X = o 

5-6X + X 2 =0 
4 



To use this program enter the following data: 

10 DATA 1, 3, 6 

11 DATA 2, 1, -6, 5 

12 DATA 3, 3, -7, 0, 4 
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NOTES: 

All 3 polynomials are solved in one run. 



Each set of coefficients must be preceded by the order of the polynomial. 

The coefficients must be entered in descending order, which is just the reverse from 
how the equations are written above. 



There is no need to put each set of coefficients on a separate line. This was done 
only for readability. 



Rooter gives erroneous results for double imaginary roots. 

SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 
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10 DATA 1. 3. 6 

11 DATA 2. 1. -6. 5 

12 DATA 3, 3, -7,~5", 4 
RUN 

WAIT. 



ROOTER 

POLYNOMIAL NUMBER 1 IS OF ORDER I 

COEFFICIENTS (IN DESCENDING ORDER) ARE: 

3 6 

THE ROOTS ARE: 
-2 

POLYNOMIAL NUMBER 2 IS OF 0RDER 2 

COEFFICIENTS [IN DESCENDING ORDER 1 ARE: 

1 -6 5 
THE ROOTS ARE: 

5 AND I 

POLYNOMIAL NUMBER 3 IS 0F ORDER 3 

COEFFICIENTS [IN DESCENDING ORDER) ARE: 

3-704 

THE ROOTS ARE: 

2 

1. AND -.666667 
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PROGRAM ABSTRACT 

PURPOSE 

To evaluate definite integrals using Gauss' rule with ten points. 

INSTRUCTIONS 

To use this program enter the functions to be integrated thus: 

10 LET Y = <A FUNCTION OF X> 

Then type "RUN. " 

The program will ask for the Integration limits as it runs. If more than one integral of the 
same function is to be evaluated, just keep supplying limits. To stop the program, make the 
two limits equal to each other. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Find the value of the following integral between the limits to 7. 65 and 7.65 to 1567. 



O log 



g 3X 



•dX 



To use this program to get a numerical approximation of the integrals type the following 
program line: 

10 LET Y = X T 2/ (LOG(3*X)/2. 30258509) 

NOTE: 

The division by 2. 30258509 was necessary to convert from natural to common logs. 
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SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

10 LET Y = XT2/(L0G(3»X>/2. 30258509) 

RW 
WAIT. 

NUMI NT 

LIMITS [L0WER, UPPER] = 7 0, 7.65 
INTEGRAL = 125.535 

LIMITS [L0WER, UPPER] : 7 7.65, 1567 
INTEGRAL = 3.64278 E 8 

LIMITS [L0VER, UPPER) = 7 0,0 
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PROGRAM ABSTRACT 



PURPOSE 

To solve spherical triangles having the apex at the North Pole and the other iwo corners 
defined by their respective latitudes and longitudes. 



INSTRUCTIONS 

Enter data for as many cases as desired successively. Data must be entered in data state- 
ments 10 - 99 in the following format: 

10 DATA LTD, LTM, LGD, LGM, RLTD, RLTM, RLGD, RLGM, ALD, ALM 

NOTE: 

Each pair of numbers specifies a location in the form Degrees, Minutes as follows: 

LTD, LTM = LOCAL LATITUDE 

LGD, LGM = LOCAL LONGITUDE 

RLTD, RLTM = REMOTE LATITUDE 

RLGD, RLGM = REMOTE LONGITUDE 

ALD, ALM = OBSERVED ALTITUDE [IF ANY] 

For South latitudes and East longitudes, enter the degree values as negative 
numbers. If there is no observed altitude, set ALD and ALM equal to zero. 

Additional instructions may be found in the listing. 



SAMPLE PROBLEM 

Solve the spherical triangle problem using the following data: 

Local Latitude - 40 deg. 50 min. North Lat. 

Local Longitude - 73 deg. 30 min. West Long. 

Remote Latitude - 23 deg. 26 min. North Lat. 

Remote Longitude - 133 deg. 30 min. West Long. 

Observed Altitude - 37 deg. 20 min. 




GE-400 SERIES TIME-SHARING 



SPHERE * * * 

M006 

Page 2 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 



10 DATA 40,50,73.30 .23.26.133.30.37.20 

rtju — — — ' — ! — l ■— 



SPHERE 

SPHERICAL TRIANGLE S L U T I N 
CASE NUMBER 1 

L0CAL P0SITI0N8 

40 DEG 50 MIN N0RTH LATITUDE 

73 DEG 30 MIN WEST L0NGITUDE 

REM0TE P0SITI0N: 

23 DEG 26 MIN N0RTH LATITUDE 

133 DEG 30 MIN WEST L0NGITUDE 

L0CAL H0UR ANGLE [AT N0RTH P0LE]: 

60 DEG 

59 DEG 60 MIN 

3 HRS 59 MIN 60 SEC 

ZENITH {GREAT CIRCLE] DISTANCES: 

52.6 DEG 

52 DEG 37 MIN 

3157 NAUTICAL MILES 

3635.5 STATUTE MILES 
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TRUE BEARINGS [GREAT CIRCLE C0URSES]: 

REM0TE P0S1TI0N FF0M L0CAL P0SITI0N: 

270.1 DEG 

270 DEG 4 MIN 

L0CAL P0SITI0N FR0M REM0TE P0SITI0N: 

55.6 DEG 

55 DEG 33 MIN 



ALTITUDE [REM0TE CELESTIAL P0SITI0N 
AB0VE L0CAL P0SITI0N H0RIZ0N): 

37.4 DEG 

37 DEG 23 MIN 



0BSERVED ALTITUDE: 

37 DEG 20 MIN 

37.33 DEG 



LINE 0F P0SITI0N: 
3 MILES AWAY 0N LINE BEARING 90.1 DEGREES TRUE 




GE-400 SERIES TIME-SHARING 



Name TRIANG * * * 
Code M007 



PROGRAM ABSTRACT 



PURPOSE 

To find the unknown features of any triangle, given one side and any two other parts. 

INSTRUCTIONS 

To use this program supply values as required by the selected option. 

NOTE: 

Specify angles as 'DEG, MIN. DEC 
Additional instructions may be found in the listing 

SAMPLE PROBLEM 

Calculate the value of the three enclosed angles of a triangle whose sides are 17, 7. 8, 
and 13.9. 

Calculate the value of three enclosed angles. 

This program has a specific way for naming the sides and angles of the triangle. Care 
must be exercised to avoid erroneous answers. 

Angle N is the angle opposite the side you choose to call side N, as shown in the following 
diagram. 




SIDE 1 
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All angles are entered as two distinct values. The first is the angle degrees and the second 
is minutes ( which can have a fractional part). If one or more of your angles is in the 
whole degrees only, you must supply for the minutes portion. 

Use the following table to select the proper option : 

If you have Us e option 

Side 1, Side 2, Side 3 1 

Side 1, Angle 3, Side 2 2 

Side 1, Side 3, Angle 1 3 

Angle 1, Side 3, Angle 2 4 

Angle 3, Angle 1, Side 3 5 



NOTE: 



The values calculated by option 1 are used to demonstrate the other options. 



SAMPLE SOLUTION 

These sample solutions are copies of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 
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THIS PR0GRAM WILL FIND THE UNKN0WN FEATURES 0F ANY 
TRIANGLE, GIVEN 0NE SIDE AND ANY TW0 0THER PARTS. 

WHAT WILL BE GIVEN tl = SSS, 2 = SAS, 3iSSA, 4=ASA, 5 = AAS]? J^ 
WHAT ARE SIDEI, SIDE2, SIDES?! 7, 7.8, 13.9 



SIDE 


17 


ANGLE [RAD] 


1 . 73268 


ANGLE tDEG) 


99 


[MINI 


16.519 


ALT T0 SIDE 


6.29426 



2 


3 


7.8 


13.9 


.469931 


.938981 


26 
55.5026 


53 

47.8778 


13.7183 


7.69802 



RADIUS 0F CIRCUMSCR CIRCLE = 8.61261 

RADIUS 0F INSCRIBED CIRCLE = 2.76492 

AREA 0F TRIANGLE = 53.5012 



AN0THER CASE [l=YES)?J^ 

WHAT WILL BE GIVEN fl =SSS, 2 = SAS, 3rSSA, 4rASA, 5 = AAS 172. 

IN0TE: SPECIFY ANGLES AS 'DEO, MIN.DEC'J 

WHAT ARE SIDEI, ANG,LE3, SIDE27 17, 53, 47.9778, 7.8 



SIDE 


17 


ANGLE [RAD] 


1.73268 


ANGLE [DEG] 


99 


[MINI 


16.5191 


ALT T0 SIDE 


6.29426 



7.8 

.469931 

26 
55.5026 

13.7183 



13.9 
.938981 

53 

47.9777 

7.69802 



RADIUS 0F CIRCUMSCR CIRCLE = 8.61261 

RADIUS 0F INSCRIBED CIRCLE = 2.76492 

AREA 0F TRIANGLE = 53.5012 



AN0THER CASE CUYESJ71 
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WHAT WILL BE GIVEN U =SSS, 2 = SAS, 3sSSA, 4= ASA, 5rAAS)?4 

[M0TEJ SPECIFY ANGLES AS *DEG, MIN.DEC) 

WHAT ARE ANG,LEI, SIDE3, AN6,LE27 99, 16.519, 13.9 , 26, 55.5026 





1 


2 


3 


SIDE 


17. 


7.8 


13.9 


ANGLE [RAD] 


1 .73268 


.469931 


.9J8981 


ANGLE [DEG] 
[MINI 


99 
16.5188 


26 
55.5024 


53 
47.9782 


ALT T0 SIDE 


6.29426 


13.7183 


7.69802 


RADIUS 0F CIRCUMSCR CIRCLE = 

RADIUS 0F INSCRIBED CIRCLE = 

AREA 0F TRIANGLE = 


8.61261 
2.76492 
53.5012 





AN0THER CASE tl=YES]?J_ 

WHAT WILL BE GIVEN [1=SSS, 2=SAS, 3=SSA, 4rASA, 5:AAS)?5 

IN0TF.: SPECIFY ANGLES AS 'DEG, MIN.DEC) 

WHAT ARE ANG.LE3, ANG.LE1, SIDE37 ^3, 47.9778, 99, 16.519, 13.9 

1 2 3 



SIDE 


17. 


ANGLE tRAD) 


1.73268 


ANGLE [DEG] 


99 


[MINI 


16.5188 


ALT T0 SIDE 


6.29426 



7.8 

.469931 

26 
55.503 

13.7183 



RADIUS 0F CIRCUMSCR CIRCLE = 8.61261 

RADIUS 0F INSCRIBED CIRCLE = 2.76492 

AREA 0F TRIANGLE = 53.5012 



13.9 

.938981 

53 
47.9776 

7.69802 
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Name CURFIT * * * 
Code M008 



PROGRAM ABSTRACT 

PURPOSE 

To determine which of 6 curves best fits the supplied data. 

INSTRUCTIONS 

Enter data on lines numbered 0-699. All values for 1 variable are given in data state- 
ments followed by all data for the second variable. One convenient arrangement might be 
lines 100-199 for variable 1 and lines 200-299 for variable 2. 

NOTE: 

This program accepts up to 200 observations on 2 variables. 

Additional instructions may be obtained by listing the program STADES * * *. 

SAMPLE PROBLEM 

Find which of 6 curves presents the closest fit to the following data points: 

1.0 3.5 3.0 7.0 

1.5 4.0 4.0 6.0 

2.0 4.5 5.0 5.0 

2.5 5.0 6.0 4.0 

3.0 5.5 7.0 3.0 

6.0 8.0 2.0 

6.5 

Supply details for curve number 3. 
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SAMPLE SOLUTION 

The following sample solution is a copy of the printout that will appear at your teletypewriter. 



NOTE: 



All user supplied information is underlined. 

100 DATA 1.1. 5. 2. 2.5. 3. 3. 5. 4. A. 5. 5. 5. 5. 6.6. 5 

200 t)ATA 3.4.5.6.7.8.7.6.5.4.3.2 
RUN 



CURFIT 



PLEASE SPECIFY THE NUMBER 0F VALUES IN] GIVEN AS DATA 
F0R THE TW0 INPUT VARIABLES, AND THE 0UTPUT C0DE ID], 
[0:1 IF 0UTPUT IS T0 BE IN 0RDER 0F INCREASING VALUES 
0F THE INDEPENDENT VARIABLE, ELSE D=0). N,D = 7 12,0 



LEAST-SQUARES CURVEFIT 
CURVE TYPE INDEX A B 



I. 
2. 
3. 
4. 

5. 
6. 



Y=A+[B*X] 

Y:A*EXPtB*X) 

Y=A*lXtBl 

Y=A+{B/X] 

Y=1/[A+B*X) 

Y:X/[A*X+B1 



5.96246 E-2 
6.79621 E-3 
4.28484 E-2 
.135872 
NEGATIVE 
NEGATIVE 



4.93421 
3.34458 
3.831 13 
2.32305 
.494583 
.328644 



E-3 



-.236842 
-4.22812 
- . 1 02 1 62 

6.09218 
-2.62454 E-2 

.148198 



F0R WHICH CURVE ARE DETAILS DESIRED [GIVE NUMBER ]?4_ 

4. Y:A+[B/X] IS A HYPERB0LIC FUNCTI0N. THE RESULTS 
0F A LEAST-SQUARES FIT 0F ITS LINEAR TRANSF0RM 
ARE AS F0LL0WS: 



X-ACTUAL 


Y-ACTUAL 


Y-CALC 


PCT DIFFER 


3 


1 


4.35378 


-77 




4 


1.5 


3.8461 


-60.9 




5 


2 


3.54149 


-43.5 




6 


2.5 


3.33842 


-25.1 




7 


3 


3.19336 


-6 




8 


3.5 


3.08457 




13.4 


7 


4 


3.19336 




25.2 


6 


4.5 


3.33842 




34.7 


5 


5 


3.54149 




41.1 


4 


5.5 


3.8461 




43. 


3 


6 


4.35378 




37.8 


2 


6.5 


5.36914 




21. 
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200 DATA 1.2,3.4,5.6,7,8,9,10.11.12 

WAIT. 



CURFIT 



PLEASE SPECIFY THE NUMBER 0F VALUES [N] GIVEN AS DATA 
F0R THE TW0 INPUT VARIABLES, AND THE 0UTPUT C0DE ID). 
t D= 1 IF 0UTPUT IS T0 BE IN 0RDER 0F INCREASING VALUES 
0F THE INDEPENDENT VARIABLE, ELSE D=0). N,D r ?I2,0 
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LEAST 
CURVE TYPE 



2. 
3. 
4. 
5. 

6. 



Y=A+[B*X) 
Y:A*EXP[B*X) 

Y = A*tXtB) 
Y=A+[B/X1 
Y=1/[A+B*X) 

Y = X/[A*X+B] 



• S U A R E S 

INDEX 

1 

.914246 

.992535 

.612308 

NEGATIVE 

.890665 



C U R V E F I T 
A B 



.5 

1.17857 

.905171 

5.1356 

.758901 

.119643 



.5 

.157214 

.771993 
-5.35806 
-6.08532 E-2 

.942427 



F0R WHICH' CURVE ARE DETAILS DESIRED [GIVE NUMBER]?! 

1. Y:A+tB*Xl IS A LINEAR FUNCTI0N. THE RESULTS 
ARE AS F0LL0WS: 



X-ACTUAL 


Y-ACTUAL 


Y-CALC 


PCT DIFFER 


1 


1 


1 


0.0 


2 


1.5 


1.5 


0.0 


3 


2 


2 


0.0 


4 


2.5 


2.5 


0.0 


5 


3 


3 


0.0 


6 


3.5 


3.5 


0.0 


7 


4 


4 


0.0 


8 


4.5 


4.5 


0.0 


9 


5 


5 


0.0 


10 


5.5 


5.5 


0.0 


11 


6 


6 


0.0 


12 


6.5 


6.5 


0.0 



F0R WHICH CURVE ARE DETAILS DESIRED [GIVE NUMBER173 

3. Y=A*[XTB) IS A P0WER FUNCTI0N. THE RESULTS 
0F A LEAST-SQUARES FIT 0F ITS LINEAR TRANSF0RM 
ARE AS F0LL0WS: 



X-ACTUAL 


Y-ACTUAL 


Y-CALC 




PCT 


DIFFER 


1 


t 


.905171 






10.4 


2 


1.5 


1.54569 




-2.9 




3 


2 


2.1138 




-5.3 




4 


2.5 


2.63947 




-5.2 




« 


3 


3.13567 




-4.3 




6 


3.5 


3.60959 




-3 




7 


4 


4.06574 




-1.6 




3 


4.5 


4.50722 




-.1 




9 


5 


4.93627 






1.2 


10 


5.5 


5.35455 






2.7 


1 1 


6 


5.76339 






4.1 


12 


6.5 


6.16382 






5.4 


<gjwh 
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Name XYPLOT * * * 
Code M009 



PROGRAM ABSTRACT 

PURPOSE 

To plot single valued functions of X, with X on the vertical axis. 

INSTRUCTIONS 

To use this program type: 

10 LET Y •= [THE FUNCTION TO BE PLOTTED] 
RUN 

During running, the program will ask for YMIN and YMAX [the limits on the horizontal 
Y-axisJ , for XMIN and XMAX [the limits on the vertical X-axis], and for DELX, the 
increment to be used along the X-axis. 

NOTE: 

Lines 11-99 of the program may be used as desired to express complicated functions. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

For illustration purposes, we will use this program to plot a normal curve. In its simplest 
form, the normal can be represented as: 

Y = e~ x ' 

From looking at the equation we can see that the Y limits are e~™ - and e° = 1. The X 
limits are ±=>, but Y remains very close to zero if X is less than -2 or greater than f2. We 
will use these as our limits and split the X axis into 40 segments ( &X - . 1). 

To use this program type in the following program line: 

10 LET Y - EX1'(-(X -A)) 
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NOTE: 



You control both the X and Y axis limits by answering the YMIN, YMAX, XMIN, XMAX, 
DELX question. 

If incorrect values were used for YMIN and/or YMAX, and the plot went off-axis at 
certain X points, the program would have typed "off-scale" and the actual Y value at 
those points. If a clean plot is required you must run the program again and modify 
YMIN and/or YMAX. 



By knowing YMIN and YMAX froni analyzing the i juation, it is possible to get a plot 
which fills all the available space. 



You can use this program to get semi-log or log- log plots by using appropriate trans- 
forms at lines 10 through 99. If this is done, though, the plot axis will still be linear 
and you must properly interpret the result. 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 
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10 LET Y = EXP(-(Xt2)) 

RUN 
WAIT. 



XYPL0T 



WHAT ARE YMIN,YMAX,XMIN,XMAX,DELX ?0,l,-2,2,.l 

F0R Xt T0P z -2 B0TT0M : 2 INCREMENT r 

F0R Y: LEFT = RIGHT = I INCREMENT : 



.1 

1.66667 E-2 



I, 



I, 



,1, 



+ 
+ 



+ 
+ 



TYPE '0* T0 ST0P 0R '1' T0 CHANGE LIMITS. WHICH70 
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Name PLOTTO * * * 
Code M010 



PROGRAM ABSTRACT 



PURPOSE 

To plot 1 to 6 functions of X simultaneously. All functions have the same upper and lower 
limits for the plot. The functions are called A, B, C, D, E, F and are plotted in that order 
of priority. Where plots would overlap, the lower priority functions are suppressed. Values 
exceeding the selected bounds are disregarded. 



INSTRUCTIONS 

To use this program enter information in the following format: 

100 DATA NUM.XMIN, XMAX, DELX, HMIN, HMAX 

200 LET A = ANY 'BASIC FUNCTION OF X 

210 LET B = ANY 'BASIC FUNCTION OF X AND/OR A 

220 LET C = DITTO FOR X AND/OR A AND/OR B 

230 [SIMILARLY FOR D] 

240 [SIMILARLY FOR E] 

250 [SIMILARLY FOR FJ 

RUN 

Where: NUM is the number of functions given [1-6], XMIN and XMAX are the 
lower and upper limits for X, DELX is the Increment for X. HMIN and HMAX 
are the lower and upper limits for the values of the function. 

NOTE: 

The horizontal increment is always HlMAX-MINj/40. 



Additional instructions may be found in the listing. 



SAMPLE PROBLEM 

Plot the same 6 curves twice where the second plot will be an enlarged version of a section 
of the first graph. 

The functions plotted are illustrated in statements 200 - 250 in the sample solutions. The 
option of plotting all 6 functions is used. Since functions D and E were out of the specified 
range, they were not shown on the second plot. 
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SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 



NOTE: 



All user supplied information is underlined. 



100 DATA 6.0,1.. 05. 0.1 

200 LET A = l-EXP(-X) ' 

210 LET B = XT2 

220 LET C = Xt3 

230 LET D = X*C 

240 LET E = EXP(-X/2) 

2 50 LET F : EXP(-X) 

RUN 

WAIT. 



PL0TT0 



MULTIPLE PL0T 0F THE FUNCTI0NS A, B, C, D, E, F 



X- VALUE 



.05 

.1 

.15 

.2 

.25 

.3 

.35 

.4 

.45 

.5 

.55 

.6 

.65 

.7 

.75 

.8 

.85 

.9 

.95 

1. 



H0RIZ INCREMENT = .025 



A 

B A 
B A 
CB A 

C B A 

DCB A 

DC B A 

DC B A 

DCB A 

DCB A 

DCB A 

DCB A 

D C B A 



FE 
F E 
F E 
F E 
F E 
F E 
F E 

F E 

F E 

F E 

F E 

F E 



DC B A F E 

D C A E 

D C F AB E 

D F C A BE 

F D A CE B 
F A D C B 

F A DC B 

F EA 



: . 





.1 



>: 
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TYPE 'S' T0 ST0P, 0R GIVE NEW VALUES 0F 
NUM,XMIN > XMAX t DELX > HMIN t HMAX?6 > 0,.5,.025 t t .l 



MULTIPLE PL0T 0F THE FUNCTI0NS A, B, C, D, E, F 







H0RIZ 


INCREMENT 


— . 


0025 



















. 


X-VALUE 



















A 














.025 


B 


A 












.05 


CB 






A 








.075 


C B 










A 




.1 


C B 












A 


.125 


DC B 














.15 


DC 


B 












.175 


D C 


B 












.2 


D C 




B 










.225 


D C 






B 








.25 


D C 








B 






.275 


D 


C 








B 




.3 


D 


C 










B 


.325 


D 




C 










.35 


D 




C 










.375 




D 




C 








.4 




D 






C 






.425 




D 








C 




.45 






D 








C 


.475 








D 








.5 










D 










• 



























.1 



TYPE 'S' T0 ST0P, 0R GIVE NEW VALUES 0F 
NUM,XMIN,XMAX,DELX,HMIN,HMAX?S 
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Name TWOPLO * * * 
Code M011 



PROGRAM ABSTRACT 

PURPOSE 

To plot simultaneously 2 functions of a single variable, X. 

INSTRUCTIONS 

To use this program type: 

10 LET Y= [THE FIRST FUNCTION OF x] 

20 LET Z= [THE SECOND FUNCTION OF X AND/OR y] 

Then: RUN 

The functions Y and Z may be any legitimate 'BASIC expressions. Intermediate variables 
may be defined using intermediate lines, if the functions are too complicated to fit on one 
line. 

Additional instructions may be obtained from the listing. 

SAMPLE PROBLEM 

Plot the first 2 curves that were used to illustrate the program PLOTTO * * *. 



NOTE: 



It is advisable to use TWOPLO * * * rather than PLOTTO * * * when only 2 functions 
are required because of the shorter running time. 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All information to be supplied by the user is underlined. 
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10 LET Y=1-EXP(-X> 
20 LET Z=Xf2 
RUN 
WAIT. 



TW0PL0 



WHAT ARE FMI N, FMAX, XHI N, XMAX, DELX70, 1 , 0, 1 , .03 



[N0TE: Y IS PL0TTED '+', Z IS '.', AND 'O' IS C0MM0N P0INT) 

F0R X: T0P : B0TT0M = 1 INCREMENT = .05 

F0R FCTS: LEFT r RIGHT = 1 INCREMENT = 1.66667 E-2 



I, 



,1. 



.1, 



,1. 



,1, 



.1. 



,1 




. + 

+ 
+ 



+ 
+ 
+ 
+ 
+ 

+ 



+ . 

+ 
+ 



TYPE 'S* T0 ST0P NSW, 0R ELSE SPECIFY NEW VALUES 
F0R FMI N,FMAX,XMIN, XMAX, DELX7S 
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Name LINPRO * * * 
Code M012 



PROGRAM ABSTRACT 



PURPOSE 



To maximize an objective function using the two-phase method. It automatically supplies 
slack, surplus, and artificial variables as required for the solution. 



INSTRUCTIONS 

NOTE: 

More detailed instructions may be obtained by listing programs LPINST * * * and 
LP1NS2 * * * 

Before using this program arrange all constraints (i. e. linear restrictions on the 
problem variables) as follows: 

1. The "less than or equal" inequalities. 

2. The strict equalities. 

3. The "greater than or equal" inequalities. 

To use this program enter the data starting at line 10,000, in this order: 

FIRST: COEFFICIENTS OF EACH OF THE PROBLEM VARIABLES (INCLUDING 

ZEROES FOR VARIABLES NOT APPEARING) IN EACH RESTRICTION 
STARTING WITH THE FIRST RESTRICTION AND PROCEEDING IN ORDER 
UNTIL ALL COEFFICIENTS OF ALL RESTRICTIONS HAVE BEEN 
ENTERED IN DATA STATEMENTS, 

THEN: THE ELEMENTS OF THE "B" VECTOR (THE CONSTANTS COMPRISING 

THE RIGHT SIDE OF ALL RESTRICTIONS) IN THE SAME ORDER AS THE 
RESTRICTIONS, 

LAST: THE COEFFICIENTS OF THE (LINEAR) OBJECTIVE FUNCTION IN THE 

SAME ORDER AS USED IN THE RESTRICTIONS, INCLUDING ZEROES 
IF NEEDED. 

THEN TYPE "RUN. " 

Additional instructions may be found in the listing. 
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SAMPLE PROBLEM 

Maximize the function 30 x XI + 45 x X2 while satisfying the following constraints: 

xi < = 6000 

X2 < = 4000 

2.5 x XI + 2.0 x X2 < = 2400 

XI > = 1000 

X2 > = 1000 

3.0 x XI - X2 > = 

2.5 x XI + 2.0 x X2 > = 10000 



The data would be entered starting at line number 10000, in the following order: 



First 

Ex 

Then 

Ex 

Last 

Ex 



Coefficients of each of the problem variables in each constraint. 

10000 Data 1, 0, 0, 1, 2.5, 2, 1, 0, 0, 1, 3, -1, 2.5, 2 
The values on the right side of the constraints. 

10001 Data 6000, 4000, 24000, 1000, 1000, 0, 10000 
The coefficients of the objective function. 

10002 Data 30, 45 



As the program runs, the values "7,2" would be supplied for "M, N. " M being the number 
of constraints and N being the number of variables. When "I, J and K" are called for 
"3, 0, 4" would be supplied. "I" being the number of < = restrictions, being the number of 
restrictions and 4 being the number of > = restrictions. 

NOTE: 

The problem solution is XI = 6000, X2 = 4000 

This problem was run 3 times to illustrate all output options. 

SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 
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10000 DATA 1.0.0.1,2.5.2,1.0.0.1.3.-1.2.5.2 

10001 DATA g50fl. f 4600. f 2J66o: lo&o! 1000^0. 10000 

10002 DATA 30. 43 

RUN 

WAIT. 



LINPR0 

TYPE '2' F0R 0UTPUT 0F TABLEAUS AND BASIS 
AT EACH ITERATION, 'I' F0R THE BASIS 0NLY, 
0R '0' F0R JUST THE S0LUTI0N. WHICH72 

WHAT ARE M AND N 0F THE DATA MATRIX?^ 



H0W ! 


fttNY 'LESS THANS* 




EQUALS' 


, 'GREATER 


THANS'73,0,4 




Y0UR 

SURPLUS 

SLACK 

ARTIFICIAL 


VARIABLES 
VARIABLES 
VARIABLES 
VARIABLES 


1 
3 
7 
10 


THR0UGH 2 
THR0UGH 6 
THR0UGH 9 
THR0UGH 13 






TABLEAU AFTER ITERATIONS 
1 
1 










6000 













1 











1 





4000 






2.5 







2 














1 
24000 






1 


0; 













-1 







1000 




1 




i 

1 




1 














-1 

1000 






3 







-1 



1 
















-1 




2.5 

-1 





2 












1 




10000 






-30 






-45 























-6.5 
1 





-2 








1 





1 





1 
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BASIS BEFORE 


ITERATION 


VARIABLE 


VALUE 


7 


6000 


8 


4000 


9 


24000 


10 


1000 


II 


1000 


12 





13 


10000 


BASIS BEF0RE 


ITERATION 


VARIABLE 


VALUE 


7 


6000 


8 


4000 


9 


24000 


10 


1000 


1 I 


1000 


1 





13 


10000 


BASIS BEF0RE 


ITERATION 


VARIABLE 


VALUE 


7 


3666.67 


8 


3000 


9 


21166.7 


10 


666.667 


2 


1000 


1 


333.333 


13 


7166.67 


BASIS BEF0RE 


ITERATION ■ 


VARIABLE 


VALUE 


7 


5000. 


8 


1000. 


9 


15500. 


4 


2000. 


2 


3000 


1 


1000 


13 


1500 


BASIS BEFORE 


ITERATION ! 


VARIABLE 


VALUE 


7 


4823.53 


8 


470.588 


9 


14000. 


4 


2529.41 


2 


3529.41 


1 


1176.47 


3 


176.471 


BASIS BEF0RE 


ITERATION ( 


VARIABLE 


VALUE 


7 


4666.67 


6 


1333.33 


9 


12666.7 


4 


3000. 


2 


4000. 


1 


1333.33 


3 


333.333 
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ANSWERS: 




VARIABLE 


VALUE 


5 


14000 


6 


13000 


9 


1000. 


4 


3000. 


2 


4000. 


1 


6000. 


3 


5000. 



DUAL VARIABLES! 

C0LUMN VALUE 

3 

4 

5 

6 

7 30. 

8 45 

9 



LDMPRO * 
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TABLEAU AFTER 6 



3 

-1 



ITERATI0NS 





2.5 







-1 




2. 
-1 





14000. 



13000. 



1 
















-2.5 







-2 






1000. 



1 


3000. 










1 











4000. 



6000. 



1 








5000. 





■I. 





30. 

8.94070 E-8 




45 

2.98023 E-8 





360000. 



4.76837 E-7 
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RUN 
WAIT. 



L1NPR0 

TYPE '2* F0R 0UTPUT 0F TABLEAUS AND BASIS 
AT EACH ITERATION, T F0R THE BASIS 0NLY, 
0R '0* F0R JUST THE S0LUTI0N. WHICH?! 

WHAT ARE M AND N 0F THE DATA MATRIX?7,Z 

H0W MANY 'LESS THANS' , 'EQUALS' , 'GREATER THANS'73,0,4 

Y0UR VARIABLES I THR0UGH 2 

SURPLUS VARIABLES 3 THR0U6H 6 

SLACK VARIABLES 7 THR0UGH 9 

ARTIFICIAL VARIABLES 10 THR0UGH 13 

BASIS BEF0RE ITERATI0N I 



VARIABLE 


VALUE 


7 


6000 


8 


4000 


9 


24000 


10 


1000 


11 


1000 


12 





13 


10000 


BASIS BEF0RE 


ITERATION 


VARIABLE 


VALUE 


7 


6000 


8 


4000 


9 


24000 


10 


1000 


11 


1000 


1 





13 


10000 


BASIS BEF0RE 


ITERATION 


VARIABLE 


VALUE 


7 


5666.67 


8 


3000 


9 


21166.7 


10 


666.667 


2 


1000 


1 


333.333 


13 


7166.67 
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BASIS BEF0RE ITERATI0N 4 

VARIABLE VALUE 

1 5000. 

8 1000. 

9 15500. 
4 2000. 

2 5000 
1 1000 
13 1500 

BASIS BEF0RE ITERATION 5 



VARIABLE 


VALUE 


7 


4823.53 


8 


470.588 


9 


14000. 


4 


2529.41 


2 


3529.41 


1 


1176.47 


3 


176.471 


BASIS BEF0RE 


ITERATION I 


VARIABLE 


VALUE 


7 


4666.67 


6 


1335.33 


9 


12666.7 


A 


3000. 


2 


4000. 


1 


1333.33 


3 


533.333 


ANSWERS: 




VARIABLE 


VALUE 


5 


14000. 


6 


13000. 


9 


1000. 


4 


3000. 


2 


4000. 


1 


6000. 


3 


5000. 


DUAL VARIABLES! 


C0LUMM 


VALUE 


3 





A 





5 





6 





7 


30. 


8 


45 


9 
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RUN 



LINPR0 

TYPE 'Z' F0R 0UTPUT 0F TABLEAUS AND BASIS 
AT EACH ITERATION, 'I' F0R THE BASIS 0NLY, 
0R '0* F0R JUST THE S0LUTI0N. WHICH70 

WHAT ARE II AND N 0F THE DATA MATRIX?7 t 2 

H0W MANY 'LESS THANS* , 'EQUALS* , 'GREATER THANS'73,0,4 



Y0UR VARIABLES 1 THR0UGH 2 

SURPLUS VARIABLES 3 THR0UGH 6 

SLACK VARIABLES 7 THROUGH 9 

ARTIFICIAL VARIABLES 10 THR0UGH 13 



ANSWERS: 




VARIABLE 


VALUE 


5 


14000. 


6 


13000. 


9 


1000. 


4 


3000. 


2 


4000. 


I 


6000. 


3 


5000. 


DUAL VARIABLES: 




COLUMN 


VALUE 


3 





4 





5 





6 





7 


30. 


8 


45 


9 
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PROGRAM ABSTRACT 

PURPOSE 

To evaluate time functions which are sums of exponentials and exponential sine- cosine terms. 

INSTRUCTIONS 

To use this program enter input data in the following format: 

20 DATA NP ,N1,N2, TO, DELTA-T, SIGMA 

where: 

NP=T0TAL NUMBER OF POINTS TO BE COMPUTED 
N1=NUMBER OF EXPONENTIAL TERMS 
N2=NUMBER OF SINE-COSINE EXPONENTIAL TERMS 
T0=TIME OF FIRST POINT 
DELTA-T=TIME BETWEEN POINTS 
SIGMA=STANDARD DEVIATION OF THE NOISE 

NOTE: 

If additive noise is not desired, SIGMA = 0. 

Parameters of the function are entered as follows: 

41 DATA c[l],c[2],c[3],... ) S[l],s[2] > s[3],... 
60 DATA a[i],a[2],a[3],...,b[i],b[2],b[3],... 
79 DATA w[l] , W[2] , w[>] , . . . , G[l] , g[Y] , G[3] , . . . 
Only statement numbers between 41 and 79 inclusive may be used. 

Maximum number of points permissable is 500 and 2N1+4N2 = 20 
The computation of Ei-X, where X is large may result in excessive running time. 
Additional instructions may be found in the listing. 
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SAMPLE PROBLEM 

Evaluate a time function of this form: 

C*EXP(-S*T) 



-0.0178189 1.5662 

0.0166119 16.2565 

0.00450666 136.889 



FORM: (A*C0S (W*T)+B*SIN (W*T) )*EXP ( -G*T) 



15.438 -.372211 491.834 200.474 



Type the following data: 

20 DATA 30, 3, 1, 0, .001, 

50 DATA -1.78189E-2, 1.66119E-2, 4.50666E-3, 1.5662, 16.2565, 136.889 

60 DATA 15.438, -.37221, 491.834, 200.474 



NOTE: 



In the data entered on line 20, 30 is the total number of points desired, 3 is the 
number of exponential terms, 1 is the number of sine-cosine exponential terms, is 
the time of the first point, .001 is the time between the points and is the standard 
deviation of the noise. 



The parameters of the functions, which are on lines 50 and 60, can be entered from 
lines 41-79. 



TMFCEV provides for adding random numbers from a uniform distribution 

(-1/2, 1/2) with a standard deviation a (standard deviation of the noise). If additive 

noise is not desired, input a as zero. 
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SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 



NOTE: 



All user supplied information is underlined. 

READY. 

0LD 

0LD PR0BLEM NAME- TMFCEV*** 

WAIT. 

READY. 

20 DATA 30.3.1.0., nni,n 
WAIT. 



TMFCEV 



I It24 



N0ISE SIGMA : 

TERMS 0F F0RM C*EXPt-S*T] ARE: 



TERMS 0F F0RM 



-1.78189 E-2 
I. 661 19 E-2 
4.50666 E-3 



1.5662 

16.2565 

136.889 



fA*C0StW*T]+B*SINtW*T])*EXPt-Q*T] ARE: 

A B w 

15.438 -.37221 491.834 



FIRST P0INT AT Tl , SPACING=T2 
TI: T2 : .001 



DATA P0INTS ARE: 

15.4413 
-4.48323 

.472038 

.329864 
-.258936 

9.37954 E-2 



10.9947 
-4.57445 

1.02218 
-.028716 
•.142646 

7.55308 E-2 



5.52152 
-3.59741 

1.30216 
-.263328 
-3.36229 E-2 

4.47195 E-2 



.603037 
-2.12686 

1.12634 
-.361716 

4.57423 E-2 

I .24639 E-2 



200.474 



•2.82 76 
-.659706 

.752094 
-.346594 

8.71895 E-2 
-.013472 



TIME: 



4 SECS. 
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PROGRAM ABSTRACT 

PURPOSE 

To provide a statistical analysis for up to 61 observations on up to 6 variables. 

INSTRUCTIONS 

To use this program enter the input data using line numbers 0-699. All values for 1 
variable are given in data statements, followed by all data for a second variable, etc. One 
convenient way to arrange the input variables is as follows: 

LINES VARIABLE 

100-199 First 

200-299 Second 

300-399 Third 

400-499 Fourth 

500-599 Fifth 

600-699 SLxth 

If no transformed variables are given a linear function is fitted for the first data variable in 
terms of the subsequent variables. By using transforms, any function of any or all of the 
input variables may be used as any variable in the analysis. Transforms are entered before 
the run, at lines 1000-1099, in the following way: 

I0XX LET X<J> = ANY FUNCTION OF V(l), V<2> , V(R) 

WHERE XCJ) IS THE J-TH VARIABLE IN THE ANALYSIS <J=1 IS THE 
DEPENDENT VARIABLE), AND THE VU) ARE THE INPUT DATA VARIABLES. 
THUS, ONE MIGHT ENTER DATA FOR THREE VARIABLES V(l>, V<2), 
AND V<3), BUT RUN THE ANALYSIS WITH X<1)=V(2> AS THE DEPENDENT 
VARIABLE, AND THREE 'INDEPENDENT' VARIABLES, BY ENTERING 
TRANSFORMATIONS LIKE THE FOLLOWING AT LINES 1000-1003: 

1000 LET X(1):V(2) 

1001 LET X(2):L0G(V(1)+V(3)) 

1002 LET X(3>=V(I)+X<2>*EXP(-1-V(3>) 

1003 LET X<4)=(V(3)-1) T 2 

Additional instructions may be obtained by listing STADES * * *. 
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SAMPLE PROBLEM 



Secure a statistical analysis for the observed data consisting of 30 data values for each of 
4 variables. 



NOTE: 

This data is identical to that used in the MULREG * * * sample problem. 

Commas are not needed after each line because this program is written in BASIC. 

The data for the first variable is entered on line 100, the second variable on 200, etc. 

The run fits the equation Y = A + B- XI + C- X2 < D-X3. The coefficients for the fitted 
equation are A = .936716, B = -.00142542, C = -.0040351 and D = . 193222. 

SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

100 DATA .43. .55. .55, .55.. 54.. 52 . .6 2. .60. .50. .55. .52, . 53. .40 

101 DATA . 46,. 57. . 59. . 49. . 45. . 48. .47. .46,. 61, .59,. 60,. 57 

102 DATA .64. .62. .54, .66, .64 

2PQ— DATA 75. 66 . 67. 70 [64. 63757.63. 7 0. 65. 73. 73. 68. 7 5. 90. 77 
201 DATA 72.70.63.66,67.57.55.51.53.58.51.51.54.54 

300 DATA 9.4.7.8.8.1 .8. 3. 8. 6. 8. 0.8 . 8. 8.7. 8.7. 8. 9. 8. 8.8. 0.8. .8. 7 

301 DATA 8.6.8.7.8.1,8.5.8.6.8.3.8.9.8.9.8.6.8.7.8.4.1.2.8.4 

302 DATA 8.3,7.8.8.5 

400 DATA 2.04.1.36.1 .36.1.36.1 .36.1.36.1.36.1 .36 t l. 59. 1.59.1 .59 

401 DATA 1.59.1.59.1.59. 1 . 1 3. 1 . 13 .1 .59 . 1 .59 . 1 .59 . 1 .59 . 1 .59. 1 .02 

402 DATA 1.02,1.02,1.13.1.13.1.13.1.36.1.13.1.13 

RUN 



MULFIT 



PLEASE SPECIFY THE NUMBER OF VALUES IN] GIVEN AS DATA 
FOR THE INPUT VARIABLES, THE NUMBER OF INPUT VARIABLES 
(R) GIVEN, THE NUMBER OF VARIABLES tS] TO BE INCLUDED 
IN THE ANALYSIS t AFTER TRANSFORMS, IF ANY], AND A CODE 
VALUE ID1. t D= 1 IF TRANSFORMS HAVE BEEN SPECIFIED, OR 
D:0 OTHERWISE.] N,R,S,D = 7 30,4,0,0 
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MULTIVARIATE C U R V E F I T 



VARIABLE REGR COEFF MEAN VALUE 

1 [CONSTANT : .936715 1 .543333 



2 
3 
4 



-J. 42542 E-3 
-4.03510 E-3 
-.193222 



65.6667 
8.24333 
1.37933 



STD DEV 

6.63995 E-2 

9.91072 
1.35663 
.241536 



STANDARD DEVIATION OF RESIDUALS = 3.43456 E-2 
INDEX OF DETERMINATION IR-SQ] : .732445 
ZERO-CHECK ON MEAN RESIDUAL = -1.24176 E-8 



ACTUAL VS CALCULATED 



ACTUAL 

.43 
.55 
.55 
.55 

.54 

.52 

.62 

.6 

.5 

.55 

.52 

.53 

.4 

.46 

.57 

.59 

.49 

.45 

.48 

.47 

.46 

.61 

.59 

.6 

.57 

.64 

.62 

.54 

.66 

.64 



CALCULATED 


DIFFERENCE 


PCT 


DIFFER 


.397707 


-3.22929 


E-2 


-8.1 




.548383 


-1.61731 


E-3 


-.2 




.545747 


-4.25326 


E-3 


-.7 




.540663 


-9.33653 


E-3 


-1.7 




.548005 


8.00545 


E-3 




1.4 


.549001 


2.90011 


E-2 




5.2 


.542922 


-7.70778 


E-2 


-14.1 




.549027 


-5.09726 


E-2 


-9.2 




.494608 


-5.39153 


E-3 


-1 




.500929 


-4.90715 


E-2 


-9.7 




.489929 


-3.00713 


E-2 


-6.1 




.493157 


-3.68432 


E-2 


-7.4 




.471776 


7.17755 


E-2 




15.2 


.487481 


2.74814 


E-2 




5.6 


.555386 


-1.46145 


E-2 


-2.6 




.573512 


-1.64875 


E-2 


-2.8 




.494179 


4.17869 


E-3 




.8 


.495415 


4.54155 


E-2 




9.1 


.50499 


2.49899 


E-2 




4.9 


.501924 


3.19242 


E-2 




6.3 


.498078 


3.80777 


E-2 




7.6 


.622468 


1 .24682 


E-2 




2. 


.62653 


3.65296 


E-2 




5.8 


.631828 


3.18277 


E-2 




5. 


.608933 


.038933 






6.3 


.630859 


-9.14134 


E-3 


-1.4 




.611784 


-8.21613 


E-3 


-1.3 




.567746 


2.77464 


E-2 




4.8 


.609929 


-5.00713 


E-2 


-8.2 




.607104 


-3.28959 


E-2 


-5.4 
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PROGRAM ABSTRACT 



PURPOSE 

To test a network of activities and events, prior to performing a critical path analysis. 
Checks are made for (1) too many activities, (2) multiple starts, (3) multiple ends, 
(4) duplicate activities, and (5) illogical loops. 



INSTRUCTIONS 

Input data is entered by activities in the following format: 

1 DATA Ij. , J x , D x , Aj , la , J a , n, , A s , ... I k , J k , D k , A* 

(Continue on lines 2-99 as needed) 

where: 

I is the starting event for the activity 

J is the ending event for the activity 

D is the duration of the activity 

A is the cost of the activity 

K is the number of activities in the network 

If no errors are detected, the program will print 

NO ERRORS. TYPE ' S' TO STOP, OR ' 1' FOR OUTPUT [~EACH 
CHARACTER IS OUTPUT INDIVIDUALLY] . WHICH? 

If a paper tape (for input to CPMTWO) is desired, the tape punch should be turned on 
immediately following the response to this question. The format of the tape is controlled 
by a data line at 130: 

130 DATA W, N, L, P 
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where: 



W is the width of the cost field (char. ) 

N is the width of the event field (char. ) 

L is the number of activities per output line 

P is a flag: P ^ requests a trace 
P = requests no trace 

As the program exists in the library, this line reads: 

130 DATA 4, 3, 3, 

so that the maximum cost for any one activity is $9999, the maximum activity number is 
999, and each output line contains 3 activities. This may be altered by the user to suit his 
particular application. 

Messages indicating errors in the network are as follows: 

1. DUPLICATE ACTIVITIES 
X I x J x D x A„ 

where: 

X is the activity number 

I x is the beginning event 

J x is the ending event 

D x is the duration of the event 

A x is the cost plus 1 . 

2. MULTIPLE STARTS 
X I x J x 

(symbols as above) 

3. MULTIPLE ENDS 
X I x J x 

(symbols as above) 
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4. LOGIC LOOP 

*1 III "x3 
^2 *X2 "x3 

Xi Ixi Jxa 

(symbols as above) 

5. TOO MANY ACTIVITIES 

RESTRICTIONS 

The starting event for the network must have only one exit, sometimes necessitating 
introduction of an extra (dummy) event. 

The program may not find all errors on a single run. Therefore, corrections should be 
made after each run, and the process continued until the NO ERRORS. . . message appears. 



SAMPLE PROBLEM 

Test the following network for errors: 




where Qj) denotes an event and » denotes an activity (Ignore cost and duration). 

SAMPLE SOLUTION 

The first run detects that the network has two end events, namely 2 and 6. Assuming that 
path from 2 to 6 has been inadvertently omitted, that path is then defined in line 12. 

The second run finds an illogical loop involving events 4, 8, and 9. Assuming that the 
direction of the activity connecting events 4 and 9 is in error, the appropriate change is 
made to line 10. 
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The third run indicates tlial the network is now error free. 

Following is a copy of the printout that will appear at your teietypewi iter. 

NOTE: 

All user supplied information is underlined. 
READY. 



1 


DATA 


3.5.0.0 


2 


DATA 


5.7.0.0 


3 


DATA 


7.2.0.0 


4 


DATA 


7.8.0.0 


5 


DATA 


8.2.0.0 


6 


DATA 


8.6.0.0 


7 


DATA 


5.9.0.0 


8 


DATA 


9.8.0.0 


9 


DATA 


8.4.0.0 


10 DATA 4*9.0*0 


1 1 DATA 4.6.0.0 



RUN 

WAIT. 

CPM0NE 



MULTIPLE ENDSl 

3 7 2 

6 8 6 

II 4 6 

12 DATA 2.6.0.0 

£UN 



CPM0NE 



L0GIC L00P» 

9 8 4 

10 4 9 

8 9 8 

9 8 4 

10 DATA 9.4.0.0 
RUN 



CPM0NE 



N0 ERRORS. TYPE ' S' T0 STOP. OR ' 1 ' FOR OUTPUT t EACH CHARACTER 
IS 0UTPUT INDIVIDUALLY!. WHICH? _S_ 
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PROGRAM ABSTRACT 



PURPOSE 

To calculate the mean, variance, and standard deviation for each of several sets of indi- 
vidual values or frequency distributions. 



INSTRUCTIONS 

To use this program enter data for each set of values as follows: 

1 DATA N, x[l], x[2j, X [3] , .x[n] 

Where the N values of the set are X [l] thru X [n] . If needed, additional data statements 
may be used to give the entire list of values. Additional cases may be given in subsequent 
data statements in the same format. 

The input for grouped values has the following format: 

1 DATA 0, N, x[l], F[l], x[2J , f[2],..., x[n] , f[n] 

Where the initial zero signals grouped data, the N is the number of different values to be 
given, and the F [i] are the number of times the X [i] occur. Data statements following 
may be used to extend the list as necessary, and blocks of grouped data may be intermixed 
freely with straight lists described above. 

As an example, suppose we were interested in the mean and standard deviation of the 
numbers 1,5,4,2,6,7,4,7 and also for the distribution consisting of 5-1'S, 3-4'S, 6-7'S, 
and 2-11 'S. These two cases could be run by typing the following: 

1 DATA 8,1,5,4,2 

2 DATA 6,7,4,7 

3 DATA 0,4,1,5,4,3 

4 DATA 7,6,11,2 
RUN 



Or eqivalently: 



1 DATA 8,1,5,4,2,6,7,4,7,0,4,1,5,4,3,7,6,11,2 
RUN 
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NOTE: 

Statement numbers 1 tliru 299 are available for continu Uion of the data field. 

Additional instruct ions may be found in the listing 

SAMPLE PROBLEM 

Determine the statistical characteristics of the following data points: 



261.4 


252.1 


255.5 


258.3 


253.2 


270.8 


268.3 


249.6 


256.3 


266.4 


265.4 


250.3 


280.9 


259.3 




261.4 


272.3 


270.3 


270.1 




258.1 


262.8 


263.2 


259.3 





To use this program to obtain the data's statistical characteristics, prepare the following 
data tape: 

1 DATA 22, 261.4, 270.8, 265.4, 261.4, 258.1, 252.1, 268.3, 250.3 

2 DATA 272.3, 262.8, 255.5, 249.6, 280.9, 270.3, 263.2, 258.3, 256.3 

3 DATA 259.3, 270.1, 259.3, 253.2, 266.4 

NOTE: 

The first data value (22) represents the total number of values. 




GE-400 SERIES TIME -SHARING 



MANDSD * * * 

S001 

Page 3 



SAMPLE SOLUTION 

These sample solutions are copies of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

READY. 

1 DAT* 22. 261.4. 270.8. 265.4. 261.4 . 258.1. 232.1 . 268,3 . 230.3 

2 DATA 272.3. 262.8 ». 255.5 . 249.6 . 280.9 , 270.3 . 263.2 . 238. I , 256.3 

3 DATA 239.3, 270.1. 259.3. 253.2, 266.4 
BJUI 

WAIT. 



MANDSD 



ARITHMETIC MEAN. VARIANCE. AND 
STANDARD DEVIATI0N 



INDIVIDUAL SET NUMBER 1 

INPUT VALUES! 261.4 
268.3 250.3 272.3 

263.2 258.3 256.3 



270.8 265.4 261.4 258.1 252.1 

262.8 255.5 249.6 280.9 270.3 

259.3 270.1 259.3 233.2 266.4 



NUMBER 0F VALUES : 22 

ARITHMETIC MEAN = 262.059 

STANDARD DEVIATI0N = 7.96691 

SAMPLE VARIANCE = 60.3866 

EST TRUE VARIANCE = 63.4716 

ST ERR0R MEAN = 1.69855 
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SAMPLE PROBLEM 



Assume that each data point has a irequency associated with it. See the following printout 
for the frequencies used. 

To use this program prepare the following data tape: 

1 DATA 0, 22, 261.4, 3, 270.8, 4, 265.4, 2, 261.4, 7, 258.1, 3 

2 DATA 252.1, > , 268.3, 7, 250.3, I, 272.3, 9, 262.8, 2, 255.5, 8 

3 DATA 249.6, 5, 280.9, 9, 270.3, 3, 263.2, 8, 258.3, 2, 256.3, 3 

4 DATA 259.3, 5, 270.1, 3, 259.3, 2, 253.2, 5, 266.4, 6 



NOTE: 

The first data value (0) indicates grouped data. The second data value (22) indicates 
the number of groups . The frequency associated with each data value is punched directly 
after that data value. 
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A g A I A 9f_ 22 ' 261 ' 4 * *' g 70 »gt j» 265.4, 2. 261.4. 7. 258.1, 3 

2"MTA 2 32.1. 3. Z6B.TTTTZS6X 1. 271.51 9. 2^8* 2, 831.8 , 8 

2 DATA 24 9.6. 5. 280.9. 9. 270.3. 3. g^j^TT T^ 8737^ 23S.~T7~T 

1DATA 23?T5, Si 270.1 , 3 t 259.3; 2; 2?3.2 |~5i~ 2«*4; € "— 

BUI 
WAIT. 



MANDSD 



ARITHMETIC MEAN, VARIANCE, AND 
STANDARD DEVIATION 



F0R GROUPED DATA SET 1 



X-VALUE 




FREQUENCY 




261.4 




3 






270.8 




4 






269.4 




2 






261.4 




7 






258.1 




3 






252.1 




5 






268.3 




7 






250.3 




t 






272.3 




9 






262.8 




2 






255.5 




8 






249.6 




5 






280.9 




9 






270.3 




3 






263.2 




8 






258.3 




2 






256.3 




3 






259.3 




5 






270.1 




3 






259.3 




2 






253.2 




5 






266.4 




6 






NUMBER 


0F VALUES 


z 


102 


ARITHMETIC MEAN 


z 


263.524 


STANDARD 


DEVIATION 


z 


8.57248 


SAMPLE 


VARIANCE 


z 


72.767 


EST 


TRUE 


VARIANCE 


z 


73.4875 




ST ERR0R MEAN 


z 


.848802 
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PROGRAM ABSTRACT 



PURPOSE 

To provide a description of univariate data with up to 300 observations on one variable. 

INSTRUCTIONS 

Enter data on lines numbered 0-699. 

Additional instructions may be obtained by listing the program STADES * * *. 

SAMPLE PROBLEM 

Determine statistical characteristics for the following data points: 



261.4 


252.1 


255.5 


258.3 


253.2 


270.8 


268.3 


249.6 


256.3 


266.4 


265.4 


250.3 


280.9 


259.3 




261.4 


272.3 


270.3 


270.1 




258.1 


262.8 


263.2 


259.3 
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SAMPLE SOLUTION 

This sample solution is a copy the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

READY. 

100 DATA 261.4,270.8,265.4,261.4,258.1.252.1,268.3.250.3.272 .3 
i~0T DATA 262. 8, 255. 5, 249. 6, 28079, 270. 3, 263. 2, 253. a/256717259. 3 
102 DATA 270.1,259.3,253.2,266.4 
RUN 



UNISTA 



PLEASE SPECIFY A TYPICAL CLASS INTERVAL 
F0R FREQUENCY DISTRIBUTIBNSi L,U =7260,280 



SUMMARY STATISTICS 

NUMBER 0F VARIATES : 22 

ARITHMETIC MEAN = 262.059 
STANDARD DEVIATI0N = 7.78374 
VARIANCE = 60.5865 
C0EFF 0F VAR IPCT] = 2.97 
STANDARD SKEWNESS = .385 
STANDARD EXCESS = -.306 



RDER STATISTICS 

SMALLEST VARIATE = 249.6 

L0WER DECILE = 250.84 

FIRST QUARTILE = 256.1 

MEDIAN : 261.4 

THIRD QUARTILE = 268.75 

UPPER DECILE = 271.85 

LARGEST VARIATE = 280.9 

T0TAL RANGE = 31.3 

DECILE RANGE = 21.01 

SEMI-QUARTILE RANGE = 6.325 

BBWLEY'S SKEWNESS = .162 

PEARS0N SKEWNESS : .254 
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FREQUENCY DISTRIBUTIBN 

UP T0 BUT PERCENT 

FROM N0T INCLUDING FREQUENCY FREQUENCY 

220 240 

240 260 10 45.455 

260 280 II 50. 

280 300 I 4.545 

300 320 



CUMULATIVE DISTRIBUTIBN 

NUMBER LESS PERCENT LESS VARIATE SUM - PCT 

VALUE THAN VALUE THAN VALUE LESS THAN VALUE 

240 

2S0 10 45.455 44.265 

280 21 95.455 95.128 

300 22 100. 100. 



0RDERED ARRAY 



249.6 


258.1 


262.8 


270.1 


250.3 


258.3 


263.2 


270.3 


252.1 


259.3 


265.4 


270.8 


253.2 


259.3 


266.4 


272.3 


255.5 


261.4 


268.3 


280.9 



256.3 261.4 
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Name ST AT AN * * * 
Code S003 



PROGRAM ABSTRACT 



PURPOSE 

To perform a statistical analysis on Data for one variable. It computes 34 difi.,rent 
measures for an array of weighted (as with frequencies) or unweighted values of the variable. 
It also gives a 10- class frequency distribution summary, and a recapitulation of the input 
data in terms of deviations from the mean and as an ordered array. 



INSTRUCTIONS 

To use this program supply data in either of these 2 formats as determined by the 
problem: 

1) FOR UNWEIGHTED VALUES: 

10 DATA 0, X(l) , X(2) , X(3) 

2) FOR DATA WITH WEIGHTS OR FREQUENCIES: 

10 DATA 1, X(l), F(l), X(2), F(2) , X(3) , F(3) 



Where the initial or 1 signals the presence or the absence of weights. Lines 11 through 
99 are available for additional input data. 



NOTES: 



This program produces output corresponding to the listing in the National Bureau of 
Standards Handbook No. 101. 



There can be no more than 99 values entered. If the values are weighted the 
maximum amount of data cannot exceed 99 values and 99 weights. 



Additional instructions may be obtained by listing the program STATEX * * *. 
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SAMPLE PROBLEM 

Suppose you had obtained the following data points as the result of an experiment: 

261.4 252.1 255.5 258.3 253.2 

270.8 268.3 249.6 256.3 266.4 

265.4 250.3 280.9 259.3 

261.4 272.3 270.3 270.1 

258.1 262.8 263.2 259.3 



To use this program to get the data's statistical characteristics you would prepare the 
following data tape: 

10 DATA 0, 261.4, 270.8, 265.4, 261.4, 258.1, 252.1, 268.3, 250.3 

11 DATA 272.3, 262.8, 255.5, 249.6, 280.9, 270.3, 263.2, 258.3, 256.3 

12 DATA 259.3, 270.1, 259.3, 253.2, 266.4 



Where the first data value (0) indicates unweighted data. 

SAMPLE SOLUTION 

These sample solutions are copies of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

10 DATA 1,261.4,3,270.8.4,265.4,2.261.4,7,258.1,3 

11 DATA 252.1.5.268.3.7.250.3.1.272.3.9.262 .8.2.255.5.8 

12 DATA 249.6.5.280.9.9.270.3.3.263.2.8.258.3.2.256.3.3 

13 DATA 259.3.5.270.1.3.259.3.2.253.2.5.266.4,6 

RUN 



WAIT. 

STATAN 

C0MPUTATI0NS 0N THE DATA ARRAY: 

NUMBER 0F VALUES = 22 

NUMBER 0F N0NZER0 WEIGHTS = 22 

SUM 0F WEIGHTS r 102 

SUM 0F UNWEIGHTED VALUES = 5765.3 

WEIGHTED MEAN : 263.524 

UNWEIGHTED MEAN = 262.059 

SMALLEST VALUE = 249.6 

LARGEST VALUE = 280.9 

RANGE : 31.3 

WEIGHTED SUM 0F SQUARES = 4.23187 E 7 

VARIANCE = 353.436 

STANDARD DEVIATION = 18.7999 

STANDARD ERR0R 0F MEAN = 1.86147 

C0EFFICIENT 0F VARIATI0N = 7.13405 E-2 
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STUDENT'S T = 141.588 

MEAN SQUARE SUCCESSIVE DIFFERENCES = 149.4 

[MEAN S9 SUCC DIFFVt VARIANCE) r .411389 

MEDIAN = 261.4 

NUMBER 0F RUNS UP AND D0WN = II 

EXPECTED NUMBER 0F RUNS = 14.3353 

STD ERROR 0F NUMBER 0F RUNS : 1.89444 

[ACTUAL RUNS - EXP RUNS1/ISTD ERR] = 1.75954 

FREQUENCY DISTRIBUTI0N [TEN EQUAL CLASSES]: 
3 2 3 .4 2 3 3 I 

COMPUTATIONS 0N DEVIATIONS FR0M MEAN: 

NUMBER 0F + SIGNS IN DEVIATIONS = 8 

NUMBER 0F - SIGNS IN DEVIATIONS = 14 

NUMBER 0F RUNS [SIGN CHANGES +11=12 

EXPECTED NUMBER OF RUNS = 11.1818 

STD ERR0R 0F NUMBER 0F RUNS = 2.10993 

(ACTUAL RUNS - EXP RUNSJ/fSTD ERR] = .387777 

TREND VALUE =-.476061 

STD ERR0R 0F TREND : 5.77038 E-2 

[TRENDl/rSTD ERR0R] =-8.25009 

BETA 0NE = 1.99451 E-4 

BETA TW0 = 8.39401 E-2 

MEAN DEVIATION = 6.58369 

RECAPITULATION OF INPUT: 



^ALUE 


DEVIATIONS 


WEIGHTS 


ORDER E; 


261 .4 


-2.12353 


3 


249.6 


270.8 


7.27647 


4 


250.3 


265.4 


1.87647 


2 


252.1 


261.4 


-2.12353 


7 


253.2 


258.1 


-5.42353 


3 


255.5 


252.1 


-11.4235 


5 


256.3 


268.3 


4.77647 


7 


258.1 


250.3 


-13.2235 


1 


258.3 


272.3 


8.77647 


9 


259.3 


262.8 


-.723527 


2 


259.3 


255.5 


-8.02353 


8 


261.4 


249.6 


-13.9235 


5 


261.4 


280.9 


17.3765 


9 


262.8 


270.3 


6.77647 


3 


263.2 


263.2 


-.323527 


8 


265.4 


258.3 


-5.22353 


2 


266.4 


256.3 


-7.22353 


3 


268.3 


259.3 


-4.22353 


5 


270.1 


270.1 


6.57647 


3 


270.3 


259.3 


-4.22353 


2 


270.8 


253.2 


-10.3235 


5 


272.3 


266.4 


2.87647 


6 


280.9 
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SAMPLE PROBLEM 



Suppose that each data point has a weight associated with it. See the printout for the 
weights used. 

To use this program prepare the following data tape: 

10 DATA 1, 261.4, 3, 270.8, 4, 265.4, 2, 261.4, 7, 258.1, 3 

11 DATA 252.1, 5, 2b8.3, 7, 250.3, 1, 272.3, 9, 262.8, :.' , 255.5, 8 

12 DATA 249.6, 5, 280.9, 9, 270.3, 3, 263.2, 8, 258. 3, 2, 256.3. 3 

13 DATA 259.3, 5, 270.1, 3, 259.3, 2, 253,2, 5, 266.4, 6 



NOTE: 

The first data value (1) indicates weighted values; the weight associated with each data 
point is punched directly after that data point. 

SAMPLE SOLUTION 

10 DATA 0, 261.4, 2 70 .8. 265.4, 261,4, 258.1. 252.1. 268.3. 230.3 

jl DATA 272.3. 262.8. 255.5. 249.6. 280.9, 270.3, 263.2, 258.3. 256.3 

12 DATA259.3. 2707T. 259.3. 253.2. 266.4 

RUN 

WAIT. 

STATAN 

C0MPUTATI0NS 9N THE DATA ARRAYt 

NUMBER 0F VALUES = 22 

NUMBER 0F N0NZER0 WEIGHTS = 22 

SUM 0F WEIGHTS = 22 

SUM 0F UNWEIGHTED VALUES = 5765.3 

WEIGHTED MEAN = 262.059 

UNWEIGHTED MEAN : 262.059 

SMALLEST VALUE = 249.6 

LARGEST VALUE = 280.9 

RANGE : 31.3 

WEIGHTED SUM 0F SQUARES = 1.51218 E 6 

VARIANCE : 63.4711 

STANDARD DEVIATI0N = 7.96688 

STANDARD ERR0R 0F MEAN = 1.69854 

C0EFFICIENT 0F VARIATI0N = 3.04011 E-2 

STUDENT'S T = 154.285 

MEAN SQUARE SUCCESSIVE DIFFERENCES : 145.4 

(MEAN SQ SUCC DI FF]/{ VARIANCE) = 2.29081 

MEDIAN : 261.4 

NUMBER 0F RUNS UP AND D0WN r 11 

EXPECTED NUMBER 0F RUNS = 14.3333 

STD ERR0R 0F NUMBER 0F RUNS = 1.89444 

[ACTUAL RUNS - EXP RUNSl/tSTD ERR] = 1.75954 




GE-400 SERIES TIME- SHARING 



STAT AN * * * 

S003 

Page 5 



FREQUENCY DISTRIBUTION [TEN EQUAL CLASSES]! 
3 2 3 4 2 3 3 1 

COMPUTATIONS ON DEVIATIONS FROM MEAN: 

NUMBER 0F + SIGNS IN DEVIATIONS = 10 
NUMBER OF - SIGNS IN DEVIATIONS = 12 
NUMBER OF RUNS {SIGN CHANGES + I ] = 12 
EXPECTED NUMBER OF RUNS = 11.9091 
STD ERROR OF NUMBER OF RUNS = 2.26883 
[ACTUAL RUNS - EXP RUNS1/ISTD ERR) = 4.00687 E-2 
TREND VALUE =-5.76505 E-2 
STD ERROR OF TREND : 6.12764 E-2 
[TREND 1/[STD ERROR J r-. 940827 
BETA ONE : .148617 
BETA TWO = 2.67928 
MEAN DEVIATION = 6.35537 

RECAPITULATION OF INPUTt 



/ALUE 


DEVIATIONS WEIGHTS ORDERE 


261. 4 


-.659089 


t 249.6 


270.8 


8.74091 


250.3 


265.4 


3.34091 


1 252.1 


261 .4 


-.659089 


1 253.2 


258.1 


-3.95909 


255.5 


252.1 


-9.95909 


256.3 


268.3 


6.24091 


258.1 


250.3 


-11.7591 


258.3 


272.3 


10.2409 


259.3 


262.8 


.740911 


1 259.3 


255.5 


-6.55909 


261.4 


249.6 


-12.4591 


261.4 


280.9 


18.8409 


262.8 


270.3 


8.24091 


263.2 


263.2 


1.14091 


265.4 


258.3 


-3.75909 


266.4 


256.3 


-5.75909 


268.3 


259.3 


-2.75909 


270.1 


270.0 


8.04091 


270.3 


259.3 


-2.75909 1 


270.8 


253.2 


-8.85909 


272.3 


266.4 


4.34091 


280.9 
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PROGRAM ABSTRACT 



PURPOSE 

To perform several types of calculations for normal and student' s-T distributions. Given 
one- or two-tailed probabilities, it will produce the corresponding limits; given a variate or 
a set of limits, it will produce the corresponding probabilities. 



INSTRUCTIONS 

To use this program enter data as follows: 

100 DATA DF, MEAN, SIGMA, (BLOCK1), (BLOCK2),. 



Where df is degrees of freedon (setting df=zero signals normal distribution), mean and 
sigma are the mean and standard deviation of the distribution, and the data blocks may be 
any mixture of the 4 types shown below: 

BLOCK(INPUT) OUTPUT PRODUCED 

1,V PROB OF X EXCEEDING THE VALUE V 

2,L,U PROB OF X INSIDE LIMITS L AND U 

3,P VALUE WHICH IS EXCEEDED WITH PROB P 

4,P LIMITS CONTAINING X WITH PROB P 



For example, the line of data shown below would produce as output the probability of a 
normally distributed variable with mean 100 and standard deviation 10 being in the range 
95-105, and also the value of the variable that will be exceeded only 2. 5 percent of the time: 

100 DATA 0,100,10,2,95,105,3,0.025 

Data for additional cases may be continued on lines 101-199 as needed. 
Additional instructions may be obtained by listing the program PROVEX * * *. 
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SAMPLE PROBLEM 

Given: Monthly sales (in pounds) will be distributed normally with a mean of 155,000. 
99 percent of the time actual sales will be within 18,000 pounds of the mean. The standard 
deviation can be approximated by 6,000 pounds. Based on the above information the following 
problems can be solved: 

1. What is the probability that sales volume will exceed 165,000 pounds? 170,000 
pounds ? 

2. What is the probability that sales volume will be at least 145,000 pounds? 

3. What is the probability that sales volume will be between 150,000 and 160,000 
pounds ? 

4. What sales volume level will fall below 25 percent of the time? 

5. What sales volume range can be expected 50 percent of the time? 

To use this program prepare the following data lines: 

100 DATA 0, 155000, 6000, 1, 165000, 1, 170000, 1, 145000, 2, 150000 

101 DATA 160000, 3, .75, 4, .5 

Question 2 is really the same as Question 1 and, therefore, uses type 1 input. 

In order to answer Question 4 above, the value which is exceeded 75 percent of the time was 
requested. 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

100 DATA 0. 155000. 6000. 1. 165000. I. 170000. I. 145000. 2. 13000 

101 DATA 160000, 3, .75, 4, .5 
RUN 



PR0VAR 



CALCULATI0NS F0R A N0RMAL DISTRIBUTION 

HAVING A MEAN 0F 155000 AND A SIGMA 0F 6000 

CASE THE PROBABILITY 0F A VARIATE: 

1 EXCEEDING 165000 IS .0478 

2 EXCEEDING 170000 IS .00621 

3 EXCEEDING 145000 IS .9522 

4 BETWEEN 150000 AND 160000 IS .59536 

5 EXCEEDING 150953. IS .75 

6 BETWEEN 150953. AND 159047. IS .5 
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SAMPLE PROBLEM 

To illustrate the use of the Student-t portion of this program, calculate several values of 
t which could be used for testing the significance between 2 sample means. Consider the 
following conditions (all of which are based on single-tail area): 

1. 20 degrees of freedom, alpha of .05 

2. 17 degrees of freedom, alpha of .025 

3. 32 degrees of freedom, alpha of .03 

To use this program type the following data lines: 

100 DATA 20, 0, 1, 3, .05 
100 DATA 17, 0, 1, 3, .025 
100 DATA 32, 0, 1, 3, .03 

Standard Student-t curves are based on a mean of and a sigma of 1. 

Since the degrees of freedom were different for each case the program was run 3 times. 

To get the two-tail areas, simply multiply by 2 the answers given by the program. 

This program will provide a critical value for 3 percent, something most tables do not provide. 

SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 
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100 DATA 20. 0. I, 3, .05 

101 

RUN 

WAIT. 



PR0VAR 

CALCULATI0NS F0R A STUDENT'S T-DISTRIBUTI0N 
HAVING A MEAN 0F AND A SIGMA 0F 1 

AND HAVING 20 DEGREES 0F FREED0M. 

CASE THE PR0BABILITY 0F A VARIATE: 

I EXCEEDING 1.72459 IS .05 



100 DATA 17. 0, 1. 3, .025 

RUN 

WAIT. 



PR0VAR 

CALCULATI0NS F0R A STUDENT'S T-DISTRIBUTI0N 
HAVING A MEAN 0F AND A SIGMA 0F I 

AND HAVING 17 DEGREES 0F FREED0M. 

CASE THE PROBABILITY 0F A VARIATE: 

1 EXCEEDING 2.10928 IS .025 



100 DATA 32, 0. 1.3. .03 

RUN 

WAIT. 



PR0VAR 

CALCULATI0NS F0R A STUDENT'S T-DISTRIBUTI0N 
HAVING A MEAN 0F AND A SIGMA 0F I 

AND HAVING 32 DEGREES 0F FREED0M. 

CASE THE PR0BABILITY BF A VARIATE: 

1 EXCEEDING 1.94981 IS .03 
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PROGRAM ABSTRACT 



PURPOSE 

To produce an equation that best fits input points described by X and Y coordinates. This 
fit is made in the least-squares sense. 



INSTRUCTIONS 

To use this linear regression program enter data using the following format: 

DATA K 

Where K represents the following functions: 

K = 1 = Linear 

K = 2 = Exponential 

K = 3 = Power 

1 DATA X[l] , Y[l] , x[2j ,y(Y] , , x[n] , y[n] 

[where x[lJ,Y[l]is the first point, x|Y] , 
Y[_2jis the second, and so on until all 
points have been entered. Additional data 
statements 2-199 may be used as needed.J 

Then type RUN. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Determine what curve best fits the following bivariate data: 

_X I X Y X Y 



1.0 
1.5 
2.0 
2.5 



38.10 
14.67 
12.76 
13.15 



3.0 

3.5 
4.0 

4.5 



11.78 
1.67 
5.35 

14.60 



5.0 

5.5 
6.0 

6.5 



5.30 

1.67 

8.92 

15.67 
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Insert the following data. This may be done by either preparing a data tape or via the 
teletypewriter. 

DATA 1 

1 DATA 1, 38.1, 1.5, 14.67, 2, 12.76, 2.5, 13.15, 3, 11.78 

2 DATA 3.5, 1.67, 4, 5.35, 4.5, 14.6, 5, 5.3, 5.5, 1.67, 6 

3 DATA 8.92, 6.5, 15.67 



NOTE: 

The data given at line indicates to fit a linear curve first. 

The exponential and power fit may be accomplished simply by entering a new line 0. 

This same data is used in the POLFIT * * * and SIXCUR * * * examples. 

SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 
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SAVE 

WAIT. 

READY. 



0LD 

0LD PROBLEM NAME- -LINREG 

WAIT. 

READY. 

DATA 1 



1 DATA 1.38. 1.1. 3.14. 67,2,12. 76.2. 5.13. IS. 3.11 . 78 

2 DATA 3.5.1.67.4^.35.4.5,n.6.5,5.3.5.T7TT6T.6 



3 DATA 8.92.6.?. 15. 67 

RUT " 

WAIT. 



LINREG 



16138 



LINEARi Y=A+B*X WITH A: 22.544 AND B=-2. 81972 



C0EFFICIEN1 


fS» C0RREL = -.527J 


239 DETERM : 


.277981 




C0MPARIS0N 


0F ACTUAL Y'S WITH 


Y'S ESTIMATED 


FR0M EQUATION: 




X-ACTUAL 


Y-ACTUAL 


Y-ESTIM 


DIFFER 


PCT-DIFF 


1 


38.1 


19.7242 


-18.3758 


-48.2304 


1.5 


14.67 


18.3144 


3.64437 


24.8423 


2 


12.76 


16.9045 


4.14451 


32.4805 


2.5 


13.15 


15.4947 


2.34465 


17.83 


3 


11.78 


14.0848 


2.30479 


19.5653 


3.5 


1.67 


12.6749 


11.0049 


658.978 


4 


5.35 


11.2651 


5.91507 


110.562 


4.5 


14.6 


9.85521 


-4.74479 


-32.4986 


5 


5.3 


8.44535 


3.14535 


59.3462 


5.5 


1.67 


7.03549 


5.36549 


321.287 


6 


8.92 


5.62563 


-3.29437 


-36.9324 


6.5 


15.67 


4.21577 


-11.4542 


-73.0966 



TIMEi 3 SECS. 
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DATA 2 



Ujs — 

WAIT. 



LINREG 



I5|3I 



EXP0NENT1A1.I Y:A*EXriB*Xl WITH A: 19.6408 AND Br-. 218326 
INDEXESt C0RREL : -.425584 DETERM = .181122 
C0MPARIS0N 0F ACTUAL Y'S WITH Y'S ESTIMATED FR0M EQUAT10N: 
X-ACTUAL Y-ACTUAL Y-ESTIM DIFFER PCT-DIFF 



1 


33.1 


15.7886 


-22.3114 


-58.5602 


1.5 


14.67 


14.1558 


•.514224 


-3.50528 


2 


12.76 


12.6918 


-6.81533 E-2 


-.534117 


2.5 


13.15 


11.3793 


-1.77069 


-13.4653 


3 


11.78 


10.2025 


-1.57749 


-13.3912 


3.5 


1.67 


9.14741 


7.47741 


447.749 


4 


5.35 


8.20143 


2.85143 


53.2977 


4.5 


14.6 


7.35327 


-7.24673 


-49.6351 


5 


5.3 


6.59283 


1.29283 


24.393 


5.5 


1.67 


5.91103 


4.24103 


253.954 


6 


8.92 


5.29973 


-3.62027 


-40.586 


6.5 


15.67 


4.75166 


•10.9183 


-69.6767 



TIME! 3 SECS* 



DATA 3 

HE 

WAIT. 



LINREG 



15x34 



P0WER1 Y=A*rXtB] WITH A= 23.5273 AMD B:-. 842416 
INDEXESt C0RREL : -.533299 DETERM = .284408 
C0MPARIS0N 0F ACTUAL Y'S WITH Y'S ESTIMATED FR0M EQUAT10NJ 



[-ACTUAL 


Y-ACTUAL 


Y-ESTIM 


DIFFER 


PCT-DIFF 


1 


38.1 


23.5273 


-14.5727 


-38.3 


1.5 


14.67 


16.7198 


2.04979 


14 


2 


12.76 


13.1214 


.361408 


2.8 


2.5 


13.15 


10.8728 


-2.27719 


-17.4 


3. 


11.78 


9.32477 


-2.45523 


-20.9 


3.5 


1.67 


8.1892 


6.5192 


390.4 


4 


5.35 


7.31792 


1.96792 


36.8 


4.5 


14.6 


6.62668 


-7.97332 


-54.7 


5. 


5.3 


6.06386 


.763862 


14.4 


5.5 


1.67 


5.59602 


3.92602 


235.1 


6. 


8.92 


5.20051 


-3.71949 


-41.8 


6.5 


15.67 


4.8614 


-10.8086 


-69.1 



TiMEl 3 SECS. 
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Namt COLINR * * * 
Code S006 



PROGRAM ABSTRACT 

PURPOSE 

To define confidence limits on simple linear regressions. 

INSTRUCTIONS 

Enter input data using the following format: 

100 DATA 116, 132, 104, 139 
200 DATA 105, 120, 85, 121 

Then type: RUN 

NOTE: 

Line numbers 0-699 are free for use as data statements, as required. All values 
for one variable are given in data statements, followed by all data for a second 
variable, and so on. One convenient way to arrange input is to put all values of 
variable 1 at lines 100-199, all for variable 2 at lines 200-299, etc. If separate 
tapes are made for each variable, the full use of input flexibility can be obtained. 

Additional instructions may be obtained by listing the program STADES * * * 

SAMPLE PROBLEM 

Determine the confidence limits using the following observed data 

Y = 116, 132, 104, 139 
X = 105, 120, 85, 121 
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SAMPLE SOLUTION 



The following sample solutions are copies of the printout that will appear at your teletype- 
writer. 



NOTE: 



All user supplied information is underlined. 



100 DATA 116. 132. 104. 139 

200 DATA 105. 120. 85. 121 

RUN 

WAIT. 



C0LINR 

SIMPLE LINEAR REGRESSI0N 
EQUATI0NJ Y = A + B*X 

H0W MANY 0BSERVATI0NS 0N EACH VARIABLE? 4 



VARIABLE 


MEAN 


VARIANCE 


STD DEVIA 


X 
Y 


107.75 
122.75 


212.688 
186.688 


14.5838 
13.6634 


INDEX fRT2] 


EXPL VAR 


UNEXPL VAR 


STD ERR0R 


.884752 


165.172 


21.5154 


4.63847 


PARAMETER 


VALUE 


95 PCT C0NF1DENCE LIMITS 


A 

B 


24.7523 
.909492 


15.8192 
.296961 


33.6853 
1.52202 



ESTIMATED VALUES 0F Y [FR0M THE REGRESSI0N) AND C0NFIDENCE 
LIMITS F0R INDIVIDUAL VALUES 0F Y, F0R EACH VALUE 0F X: 



X-ACTUAL 

105 
120 
85 
121 



Y-ACTUAL 

116 
132 
104 
139 



Y-CALC 

120.249 
133.891 
102.059 
134.801 



95 PCT C0NFIDENCE LIMITS 



100.203 
112.554 
77.7037 
1 13.24 



140.295 
155.229 
126.414 
156.362 



TIMEj 3 SECS. 
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SAMPLE PROBLEM 

Determine the confidence limits using the following observed data: 

Y = -3, -2, -1, -0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 
X = 10, 9, 8, 7, 6, 5, 4, 3, 2, 1, 0, -1, -2, -3 



SAMPLE SOLUTION 

100 DATA 'I, -2, -I. -0, I, 2, 3, 4, 5, 6, 7, g, 9, 10 
200 DATA 10, 9, 8, 7, ^ 7~ 5» 4 » 3 » 2 » '» °t '!■ -2. -3~ 

RTJB 

WAIT. 



C0LINR 

SIMPLE LINEAR REGRESSION 
EQUATION] Y = A + B*X 

H0W MANY OBSERVATIONS 0N EACH VARIABLE? 14 
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VARIABLE 

X 
Y 



MEAN 

5.5 
3.5 



VARIANCE 

16.25 
16.25 



STD DEVIATION 

4.031 13 
4.031 13 



INOEX (RT2) EXPL VAR UNEXPL VAR STD ERROR 

1 16.25 



PARAMETER 

A 
B 



VALUE 

7 
•I 



95 PCT CONFIDENCE LIMITS 



7 
-I 



ESTIMATED VALUES OF Y [FROM THE REGRESSION) AND CONFIDENCE 
LIMITS FOR INDIVIDUAL VALUES OF Y, FOR EACH VALUE OF X: 



X-ACTUAL 


Y-ACTUAL 


Y-CALC 


10 


-3 


-3 


9 


-2 


-2 


8 


-1 


-I 


7 








6 


1 


1 


5 


2 


2 


4 


3 


3 


3 


4 


4 


2 


5 


5 


1 


6 


6 





7 


7 


-1 


8 


8 


-2 


9 


9 


-3 


10 


10 


TIME: 4 


SECS. 




"($fmk, 






Wuf&f 







95 PCT CONFIDENCE LIMITS 



-3 
-2 
-1 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 



-3 
-2 
-t 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 



GE-400 SERIES TIME- SHARING 



Name CONBIN * * * 
Code S007 



PROGRAM ABSTRACT 



PURPOSE 

To determine confidence limits for a population proportion based on the normai curve 
approximation. 

NOTE: 

For small sample sizes use the program BICONF * * *. 



INSTRUCTIONS 

To use this program enter values for X and N when requested by the program. 



NOTE: 



X = number of successes in sample 
N = sample size 



Additional instructions may be found in the listing. 



SAMPLE PROBLEM 

A polling agency makes a sample of 200 voters in a certain city and it is found that 110 of 
these people intend to vote for Candidate A. Suppose the agency is willing to settle for 
accuracy of only 50 percent (or 75, or 90, or 95, etc. ) in their prediction, what limits 
should they choose? 

Suppose that they decide to make a larger survey for greater accuracy. This time they poll 
2000 people and find that 1150 people (or 57. 5 percent) are planning to vote for Candidate A. 
What limits should they then pick for a 50 percent (or 75, or 90, or 95, etc. ) chance of 
accuracy? 
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SAMPLE SOLUTION 

These sample solutions are copies of the printout that will appear at your teletypewriter. 



NOTE: 



All user supplied information is underlined. 



RUN 
WATT. 



CONBIN 



CONFIDENCE LIMITS FOR A POPULATION PROPORTION BASED ON 
THE NORMAL CURVE APPROXIMATION. IFOR SMALL SAMPLE SIZES 
USE THE PROGRAM BICONF***!. WHAT ARE X [NUMBER OF SUCC- 
ESSES IN SAMPLE), N tSAMPLE SIZE17 IIP, 200 

BEST ESTIMATE OF POPULATION PROPORTION IPCT1 = 55. 



CONFIDENCE LIMITS ON POPULATION PROPORTION 
CONF LEVEL LOWER LIM UPPER LIM 



50 


52.377 


57.623 


75 


50.703 


59.29 7 


90. 


48.964 


61.036 


95. 


47.855 


62.145 


99. 


45.689 


64.311 


99.9 


43.174 


66.826 


99.99 


41.063 


68.937 


RUN 






WAIT. 







CONBIN 

CONFIDENCE LIMITS FOR A POPULATION PROPORTION BASED ON 
THE NORMAL CURVE APPROXIMATION. IFOR SMALL SAMPLE SIZES 
USE THE PROGRAM BICONF***!. WHAT ARE X I NUMBER OF SUCC- 
ESSES IN SAMPLE!, N (SAMPLE SIZE!? 1150, 2000 

BEST ESTIMATE OF POPULATION PROPORTION tPCT) = 57.5 

CONFIDENCE LIMITS ON POPULATION PROPORTION! 
CONF LEVEL LOWER LIM UPPER LIM 



50 


56.729 


58.271 


75 


56.203 


58.797 


90. 


55.657 


59.343 


95. 


55.308 


59.692 


99. 


54.628 


60.372 


99.9 


53.838 


61.162 


99.99 


53.174 


61.826 
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Name CONDIF * * * 
Code S008 



PROGRAM ABSTRACT 



PURPOSE 

To compute confidence limits for the difference between 2 population means, : .,ased on data 
supplied for 2 samples, 1 from each population. 



INSTRUCTIONS 

To use this program type the data in the following format: 

100 DATA HI, Nl, Ml, SI, H2,N2,M2,S2 
THEN RUN 

NOTE: 

HI = SIZE OF POPULATION 1 [LET HI 

EQUAL ZERO IF POPULATION IS INFINITE] 

Nl = SIZE OF SAMPLE 1 

Ml = ARITHMETIC MEAN OF SAMPLE 1 

SI = STANDARD DEVIATION OF SAMPLE 1 
[BASED ON DIVISOR OF Nl] 

H2, N2, M2, S2 ARE THE SAME FOR SAMPLE 2. 
Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Determine the confidence limits for these 2 samples 

(1) 0, 21, 30.9048, 29.2752 

(2) 0, 21, 44, 26.8239 




GE-400 SERIES TIME- SHARING 



CONDIF * * * 

S008 

Page 2 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 



NOTE: 



All user supplied information is underlined. 



100 DATA 0.21. 30. 9048. 29.2732. 0.21. 44. 26. 8239 

RUN 

WAIT. 



C0NDIF 



STATISTIC 



SAMPLE 1 



SAMPLE 2 



SAMPLE MEAN 


30.9048 




44 


SAMPLE VARIANCE 


857.037 




719.522 


SAMPLE STD DEVIATI0N 


29.2752 




26.8239 


SAMPLE SIZE 


21 




21 


P0PULATI0N SIE 


INFINITE 




INFINITE 


ESTIM P0PN STD DEV 


29.9982 




27.4863 


STD ERR0R 0F MEAN 


6.54613 




5.99801 


DIFF BETWEEN MEANS 




-13 


.0952 


STD ERR0R 0F DIFF 




8.1 


37851 


DEGR 0F FREED0M IDIFF] 




39 


.7 



C0NFIDENCE LIMITS 0N DIFFERENCE BETWEEN MEANS: 



C0NF LEVEL 


L0WER LIM 


UPPER LI 


50 


-19.1389 


-7.05145 


75 


-23.4602 


-2.73022 


90. 


-28.0478 


1.85743 


95. 


-31.0433 


4.85287 


99. 


-37.1143 


10.9239 


99.9 


-44.6336 


18.4432 


99.99 


-51.4662 


25.2758 


99.999 


-57.9524 


31.762 
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Name CONLIM * * * 
Code S009 



PROGRAM ABSTRACT 



PURPOSE 

To compute confidence limits for an unknown population mean based on random sample 
data. 



INSTRUCTIONS 

To use this program enter the input data using this format: 

50 DATA [SIZE OF POPULATION] 

[omit THIS INPUT IF INFINITE pop'n] 
100 DATA xQlJ,X^2j ,X Qn ] 

WHERE THE X [l ] ARE THE SAMPLE OBSERVATIONS. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Analyze the data which is entered in the following sample solution. 




GE-400 SERIES TIME-SHARING 



CONLIM * * * 

S009 

Page 2 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 



NOTE: 



All user supplied information is underlined. 



100 DATA 2. 3. 6. 9.8. It. 13. 16. 23. 45. 67 .89.35. 67. 13. 43. 21. 67. 98. 7. 6 

5Tm 



CONLIM 



VALUES OF SAMPLE STATISTICS » 

SIZE OF SAMPLE 21 

SAMPLE MEAN VALUE 30.9048 

VARIANCE OF SAMPLE 857.039 

SAMPLE STD DEVIATION 29.2752 

ESTIMATED POPN STD DEV 29.9982 

STANDARD ERROR OF MEAN 6.54614 



CONFIDENCE LIMITS ON POPULATION MEANt 
CONF LEVEL LOWER LIM UPPER LIM 



50 


26.4895 


35.32 


75 


23.3745 


38.4351 


90. 


20.1373 


41.6722 


95. 


18.0745 


43.735 


99. 


14.0428 


47.7667 


99.9 


9.36435 


52.4452 


99.99 


5.43601 


56.3735 


99 .999 


1 .99002 


59.8195 
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Name BICONF * * * 
Code S010 



PROGRAM ABSTRACT 



PURPOSE 

To determine the confidence limits for a population proportion based on the e. ct binomial 
distribution. 



INSTRUCTIONS 

Enter values for X, N, and C when requested by the program. 

NOTE: 

X - Successes. 

N - Sample size. 

C - Confidence Coefficient in percent. 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

A polling agency makes a sample of 200 voters in a certain city and it is found that 110 of 
these people intend to vote for Candidate A. Therefore, the best estimate that can be made 
from this sample is that 55 percent of the entire population intend to vote for Candidate 
A. 

If the agency wants to publish a prediction, with a 95 percent chance that they will be 
correct that the actual percentage of the entire population will be within certain bounds, 
what limits should they choose? 
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SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

RUN 
WAIT. 

BIC0NF 

CONFIDENCE LIMITS FOR A POPULATION PROPORTION BASED ON 
THE EXACT BINOMIAL DISTRIBUTION. WHAT ARE THE VALUES 
0F XtSUCCESSES], NtSAMPLE SIZE], CtCONFIDENCE COEFF- 
ICIENT IN PERCENT)? ! 10,200,95 

BEST ESTIMATE OF POPULATION PROPORTION [PCT) = 55. 

THE 95 PERCENT CONFIDENCE LIMITS ON THE POPULATION 
PROPORTION [PCT] ARE 47.8249 AND 62.0247 
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Name BITEST * * * 
Code SOU 



PROGRAM ABSTRACT 

PURPOSE 

To perform a statisical test of a binomial proportion. 

INSTRUCTIONS 

Enter values for X, N, and P when requested by the program. 

NOTE: 

X - Successes in sample 

N - Sample size 

P - Population proportions 

Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Consider a city in which 75 percent of the population intend to vote for Candidate A (and the 
rest for some other candidate). From a survey of 200 people picked at random, what is the 
probability that 60 percent or less (i.e. 120 people) are planning to vote for Candidate A? 

Let a "success" be a person in the sample who intends to vote for Candidate A. Therefore, 
the input to the program will be: 

X(Number of successes in sample) = 120 
N( Sample size) = 200 

P(The true proportion of the 

population intending to vote 

for A) = . 75 

As can be imagined, the accuracy of a smaller sample (say 20 people instead of 200) is much 
less. This is demonstrated by the second of the 2 sample solutions. 
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SAMPLE SOLUTION 

These sample solutions are copies of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

READY. 

RUN 
WAIT. 

BITEST 

THIS PR0GRAM MAKES THE NECESSARY CALCULATI0N F0R A 
STATISTICAL TEST 0F A BIN0MIAL PR0P0RTI0N. WHAT 
ARE (SUCCESSES IN SAMPLE), NtSAMPLE SIE], PtTHE 
P0PULATI0N PR0P0RTI0N)? I2O,2QO,.75 

IN SAMPLES 0F SIE 200 RAND0MLY SELECTED FR0M A 
B1N0MIAL P0PULATI0N HAVING A TRUE PR0P0RTI0N 0F .75 
THE PR0BABILITY 0F A SAMPLE HAVING 120 0R LESS 
SUCCESSES IS .000002 



RUN 
WAIT. 



BITEST 



THIS PR0GRAM MAKES THE NECESSARY CALCULATI0N F0R A 
STATISTICAL TEST 0F A BIN0MIAL PR0P0RTI0N. WHAT 
ARE (SUCCESSES IN SAMPLE], NtSAMPLE SIE), PITHE 
P0PULATI0N PR0P0RT10N]? 12,2Q,.75 

IN SAMPLES 0F SIE 20 RAND0MLY SELECTED FR0M A 
BIN0MIAL P0PULATI0N HAVING A TRUE PR0P0RTI0N 0F .75 
THE PR0BABILITY 0F A SAMPLE HAVING 12 0R LESS 
SUCCESSES IS .101812 
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Name BINDIS * * * 
Code S012 



PROGRAM ABSTRACT 



PURPOSE 

To compute the binomial distribution function for problems of this type: In a series of N 
independent trials, each with probability of success P, what is the probability that there 
will be exactly X successess? 



INSTRUCTIONS 

Enter values for P,N, and M, when requested by the program. 

NOTE: 

This program will determine all the values of X for which the probability is sign- 
ificantly large [up to 7 standard deviations away from P * N] 

The number M is a multiplier of the Binomial function useful in scaling the function to a 
conveniently graphable range. [M does not affect the values of the binomial function, but 
is multiplied and printed out in a separate column. Each value of B *M is rounded off to 
the nearest integer] 

Additional instructions may be found in the listing. 



SAMPLE PROBLEM 

Determine the probability that in 30 rolls of a die it will come up 3 exactly 6 times (or 7, or 
8, or 9, etc). 



To use this program enter the following data: P = 1/6 = . 166667 

N =30 
M = 100 
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SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 



NOTE: 



All user supplied information is underlined. 



RUN 
WAIT. 



BINDIS 

D0 Y0U WANT INSTRUCTI0NS (I = YES, O=N01? 

WHAT ARE P,N,M? .166667, 30, 100 



2 

3 

4 

5 

6 

7 

3 

9 

10 

11 

12 

13 

14 

15 

16 

17 



IUTI0N: 


N : 


30 P : 


.166667 


X/N 




BtX,P,Nl 




B * 







4.21267 


E-3 





3.33333 


E-2 


2.52761 


E-2 


3 


6.66667 


E-2 


7.33008 


E-2 


7 


.1 




.136828 




14 


.133333 




.184719 




18 


.166667 




.192108 




19 


.2 




.16009 




16 


.233333 




.109777 




II 


.266667 




6.31217 


E-2 


6 


.3 




3.08596 


E-2 


3 


.333333 




1.29611 


E-2 


1 


.366667 




4.71312 


E-3 





.4 




1.49249 


E-3 





.433333 




4.03306 


E-4 





.466667 




1.00375 


E-4 





.5 




2.14133 


E-5 





.533333 




4.01501 


E-6 





.566667 




6.61297 


E-7 






100 
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Assembly Total Units Number Percent 

Line Assembled Rejected Rejected 



1 1217 45 3.7 

2 948 49 5.2 

3 1165 33 2.8 

4 1121 44 3.9 



Line 2 seemed higher, and line 3 lowest in reject rate, but such a difference could be the 
result of just chance. The engineer decided to make a statistical test for significance. 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 

READY. 

100 DATA 1217, 45, 948, 49. 1165, 33. 1121. 44 

RUN" 

WAIT. 

SEVPR0 11:14 

CHI-SQUARE TEST 0F SEVERAL PR0P0RTI0NS: 
SAMPLE SUCCESSES / T0TAL PCT SUCCESSES 



1 


45 


/ 1217 


3.698 


2 


49 


/ 948 


5.169 


3 


33 


/ 1 165 


2.833 


4 


44 


/ 1121 


3.925 



CHI-SGUARED F0R THESE DATA = 7.81993 
C0RRESP0NDING N0RMAL DEVIATE : 1.71866 

THE PR0BABILITY 0F THIS VALUE 0F CHI-SQUARE 
BEING EXCEEDED BY CHANCE AL0NE IS .043 



TIME: 2 SECS. 



ANALYSIS OF RESULT 

It looks as though the 4 lines are not alike on quality. Lines 1 and 4 are much alike, but 
Line 2 seems to have problems. The test of the four proportions simply tells us that it's 
rather improbable (only 1 chance in 24 or so) that chance could account for this much 
variability, assuming the quality level was actually the same on all lines. 
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Name SEVPRO * * * 
Code S013 



PROGRAM ABSTRACT 

PURPOSE 

To apply a CHI-SQUARE test to several sample proportions. 

INSTRUCTIONS 

To use this program simply supply values for N and S in the following format: 



100 DATA N1.S1, N2.S2, N3,S3, 
RUN 



where: 

Nl = size of sample 

SI - the number of successes in sample 1. 

Repeat for the total number of samples. 
Additional instructions may be found in the listing. 

SAMPLE PROBLEM 

Statistically test the following data for significance. 

In response to heavy demand for a particular model portable radio, four separate assembly 
lines have been in operation for the last two weeks. While all are identical operations for 
all intents and purposes, there are unavoidable differences in equipment, operator experience 
and so on. The reject rate has been running fairly high, and each line is blaming it on the 
others. The Quality Engineer decided to check all rejects for one day to discover whether 
the quality was significantly different for the four lines. This data appears on the following 
page. 
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Name ONEWAY * * * 
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PROGRAM ABSTRACT 



PURPOSE 

To perform a one way analysis of variance with equal sample sizes. 

INSTRUCTIONS 

To use this program enter data using the following format: 

100 DATA M,N 

200 DATA X11,X12,X13, ...,X1N 

201 DATA X21,X22,X23,...,X2N 
ETC UNTIL 

2XX DATA XM1,XM2,XM3,.. . ,XMN 

where: 



M = number of samples 
N ~ common sample size 
XIJ = observations 



Additional instructions maybe found in the listing. 



SAMPLE PROBLEM 

Determine what the probability is that the variability in the observed data is due to pure 
chance. 

As Manager of the Materials Laboratory, you have been following with interest the 
preparation of tensile test specimens on three rather exotic alloys being developed for 
high-temperature use. The materials are difficult to fabricate, but the men have finally 
succeeded in getting four test bars ready for stretching in the tensile test machines. Since 
these alloys differ only in the relative amounts of a single ingredient, the main purpose for 
the test is to discover the effect of the ingredient on the tensile strength, if indeed it has 
any at all. Results of the twelve tests are on the following page. 
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(Tensile Values in 1000 lbs./sq. inch) 



Alloy A 


Alloy B 


Alloy 


181.4 


166.3 


176.4 


168.6 


172.1 


183.2 


174.0 


179.2 


189.4 


183.2 


159.7 


177.5 



What's the probability that the observed variability is due to pure chance rather than the 
amount of the added ingredient? 



SAMPLE SOLUTION 

This sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined. 
100 DATA 3^ 4 



200 DATA 181.4, 168.6, 174, 183.2 

201 DATA 166.3. 172.1. 179.2. 159.7 

202 DATA 176.4. 183.2. 189.4. 177.5 

RUN 
WATT. 



0NFWAY 11:30 



ANALYSIS F VARIANCE 
3 SAMPLES 0F SIZE 4 

SAMPLE SAMPLE T0TAL SAMPLE MEAN 

1 707.2 176.8 

2 677.3 169.325 

3 726.5 181.625 

MEAN SQUARE [BETWEEN SAMPLES) = 153.631 

MEAN SQUARE [WITHIN SAMPLES J : 50.1616 

CALCULATED VALUE 0F F-RATI0 : 3.06273 

C0RRESP0NDING N0RMAL DEVIATE : 1.30578 

THE PR0BABILITY 0F AN F-RATI0 THIS LARGE 
0CCURRING BY CHANCE AL0NE IS .096 



TIME: 2 SECS. 



ANALYSIS OF RESULT 

The evidence as to the effect of varying amounts of the alloying ingredient is not too clear. 
If the ingredient had no effect at all, and we ran ten experiments like this, we'd expect one 
of the ten to show as much diversity of sample means as we actually found. Maybe the test 
method should be looked at for improvement; there's quite a bit of variability among 
specimens of the same alloy. 
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Name TESTUD * * * 
Code SOI 5 



PROGRAM ABSTRACT 

PURPOSE 

To test an unknown population mean using sample statistics. 

INSTRUCTIONS 

To use this program simply supply values for the 5 variables N, M, S, W, andX. 

where: 

N = Sample size 

M = Sample mean 

S = Sample standard deviation 

W = Population size (0 if infinite) 

X = The population mean to be tested 



Additional instructions may be found in the listing. 



SAMPLE PROBLEM 

Determine the probability that the annual fallout will exceed or equal a critical amount. 

During the early stages of a project to develop high-energy fuels for gas turbine engines, the 
question of liability for damages to agricultural crops in the vicinity of the outdoor test site 
arose. As the exhaust products were toxic to plant growth if applied in concentrations 
exceeding about 20 pounds per acre, a sampling experiment was established to measure 
the fallout at various distances from the test site. The greatest potential liability seemed 
to be from farmland located about a mile and a half downwind from the prevailing direction 
of local winds, since crop losses for any reason would doubtless be blamed on the "poison 
gases" which the local people were already grumbling about. By annualizing the results 
from the first eight test runs, a sample of eight readings averaging 13.4 pounds per acre 
with a standard deviation of 5.1 p.p.a. was available for the lawyers' consideration. Their 
question was simply this: is this evidence sufficient to deny a claim that the fallout actually 
will equal or exceed the critical value of p.p.a. ? 
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SAMPLE SOLUTION 

'Phis sample solution is a copy of the printout that will appear at your teletypewriter. 

NOTE: 

All user supplied information is underlined 
RUN 



TESTUD 

THIS PR0GRAM PERF0RMS CALCULATI0NS NECESSARY F0R 
TESTING AN UNKN0WN P0PULATI0N MEAN USING SAMPLE 
STATISTICS. WHAT ARE N {THE SAMPLE SIZE), M [THE 
SAMPLE MEAN], S [THE SAMPLE STANDARD DEVIATI0N), 
W [P0PULATI0N SIZE, ZEF0 IF INFINITE], AND X [THE 
P0PULATI0N MEAN T0 BE TESTED17 8, 15.4, 5.1, 0, 20 

BASED 0N THE STUDENT'S T-DISTRIBUTI0N WITH 7 
DEGREES 0F FREED0M, THE PR0BABILITY 0F FINDING 
A SAMPLE MEAN THIS MUCHLESS THAN THE P0PUL- 
ATI0N MEAN IS .0006 

TIME: 2 SECS. 

ANALYSIS OF RESULT 

In terms of the statistical model, we conclude that such sample results as we observed would 
be extremely rare if the population mean were 20 p.p. a. It seems fairly safe to say that the 
annual fallout will be less than this critical amount. 




GE-400 SERIES TIME-SHARING 




UTHQ U.8A 



